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Knowledge  ol  whether  or  not  all  modes  ol  consequence  have  been  found  is  ol  paramount  importance  in 
numerical  waveguide  studies  A general  loot  tinding  method  (Rl  M).  developed  Tor  vlf/lt  propagation  in  the  earth- 
ionosphere  waveguide,  locates  all  complex  zeros  trt  physically  important  rectangular  regions  ol  lire  complex  plane. 
The  method,  based  on  phase  contour  tracing,  is  applied  to  the  problem  ol  tropospheric  ducting  in  a trihncai 
reluctivity  environment  frequencies  Irom  vlit  and  into  the  microwave  range  are  considered.  Principles  upon 
which  the  RI  M is  founded  are  reviewed  and  modifications  ol  the  modal  equation  demanded  by  the  method  are 
discussed  Particularly  in  the  microwave  range,  overflow  can  he  a problem  and  methods  of  overcoming  this  are 
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discussed  llowevet.  ( Ik'  KI  M is  most  easily  implemented  il  the  latter  modifications  can  he  avoided  Tlierelore. 
tables  aie  included  to  indicate  lot  trilmeai  pioliles  the  maximum  liecpiency  below  which  no  oveillow  problems 
would  occui  riiese  tables  depend  upon  the  maximum  number  which  the  computer  can  handle,  the  location  id 
the  search  lectangle.  and  the  ref ractivity  profile  parameters  In  the  present  case  stud)  ot  duct  nip  produced  by  a 
stione  elevated  layer  (40  M unit  delkit  at  about  ’()()  nr)  the  KI  M has  been  applied  successfully  into  S band  with 
computer  run  tune  reejini emei 1 1 s computable  to  those  ol  the  more  conventional  Newton  Kaphsou  method 
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OBJECTIVE 


I \plore  the  u 1 1 1 1 1 \ ot  .1  recently  developed  roul  finding  method  ( R I M ) lor  locating 
waveguide  modal  solutions  associated  with  tropospheric  ducting  produced  In  trilinear  re- 
tractiv  ity  env iromnents 


RESULTS 

In  the  present  case  study  ot  ducting  produced  by  a strong  elevated  layer  (40  M unit 
deficit  at  about  200  m)  the  RIM  has  been  applied  successfully  from  vhl  into  S band  with 
computer  run  time  requirements  comparable  to  those  ot  the  more  conventional  Newton- 
Raphson  method.  I he  RIM.  however,  unlike  the  Newlon-Raphson  method,  assures  the 
location  ot  all  significant  modes 


RECOMM  ENDATION 

Investigate  the  utility  ot  the  RIM  lor  locating  modes  in  multisegmented  tropospheric 
retractiv  ity  environments 
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I INTRODUCTION 


It  is  well  known  that  the  presence  of  strong  gradient,  layered  structure  ot  the  tropo 
sphere  greatly  influences  the  propagation  of  radio  waves  ol  frequency  greater  than  30  Mil/ 

(and  perhaps  lower  under  strong  layering  conditions)  to  the  extent  that  beyond  the  horizon 
ducted  fields  tnav  be  main  tens  ol  dB  above  the  troposcattci  lields  It  has  been  recognized 
lor  a long  time  that  waveguide  concepts  can  be  used  to  at  least  qualitatively  explain  much 
anomalous  propagation  and  with  the  advent  ol  the  modern  day  computer  a number  ol  rather 
recent  studies''*’  have  been  concerned  with  rather  detailed  calculations  in  tropospheric  duct- 
ing environments  based  on  waveguide  concepts.  Waveguide  theory  treats  the  held  as  being 
composed  of  one  or  more  discrete  families  (modes)  of  plane  waves  conlined  to  a waveguide 
which,  in  the  case  ot  tropospheru  ducts,  has  an  upper  boundary  that  can  vary  between  tens 
of  metres  and  thousands  of  metres  The  number  of  waveguide  modes  needed  to  adequately 
represent  the  field  for  tropospheric  ducts  produced  by  elevated  layers  can  typically  range 
Irom  one  mode  at  the  lower  end  of  the  v hi  band  to  hundreds  ot  modes  at  S band  I he 
propagation  parameters  tie.  phase  velocities  and  attenuation  rates)  depend  upon  the  geometry 
ol  the  guide  and  principally  upon  the  tropospheric  rel ractivity  (or  refractive  index)  profile. 

I or  a horizontally  stratified  guide,  once  all  significant  modes  have  been  found,  the  radio 
lields  can  be  determined  for  a given  transmitter-receiver  configuration  by  simply  summing 
(with  proper  weights)  the  significant  modes. 

Unfortunately,  for  realistic  layered  refractivity  profiles,  solutions  to  the  modal  equa- 
tion cannot  be  found  explicitly  but  rather  must  be  found  by  iterative  techniques  so  that  a 
crucial  problem  in  numerical  waveguide  studies  is  knowing  whether  or  not  all  modes  ol 
consequence  have  been  found.  .Mathematically  this  reduces  to  the  problem  ot  finding  the 
complex  zeros  of  a modal  equation  in  some  region  of  a complex  eigenvalue  space.  A 
quotation  from  Hamming  concerning  the  problem  ol  locating  complex  zeros  ol  a function 
seems  appropriate  at  this  point. 

"It  is  curious  that  the  real-zero  problem  has  been  extensively  investigated  while  the 
complex-zero  problem  has  generally  been  ignored  I vidently  it  is  a field  ripe  lor 
lurther  research." 

One  recent  root  finding  method^  ( RI  M I developed  primarily  for  vll  and  lower  11 
( 10  kHz  to  about  00  kHz)  propagation  in  the  earth-ionosphere  waveguide  finds  all  modes  in 

1 Wait.  JR.  and  Spies,  kl’.  Internal  guiding  ol  microwaves  by  an  elevated  tropospheric  layer.  Radio  Sci 
vol  4 l‘)(>x».  p .H9-d2b 

2 ( hang.  Ill . I he  etleet  ol  tropospheric  layer  structures  on  long  range  VIII  radio  propagation.  II  hi  Trans 

Antennas  I'ropagal . API9,  751-756 

4  Dresp.  MR.  Radio  wave  propagation  in  the  presence  ot  an  elevated  tropospheric  duel,  PhD  diesis  ( 1972) 
t mversit)  ol  Pennsylvania.  Philadelphia.  Pennsylvania 

4 Dresp  MR.  I ron1 2 * 4 5 6 7  spheric  duct  propagation  al  VIII  . t ill  and  Sill  . Mitre  Corporation  I echnical  Report 
MIR  41  14.  vol  I.  October  1975 

5 Papperi  KA.  and  Goodhart.  CL.  ( ase  studies  ol  beyond-the-hori/on  propagation  in  tropospheric  ducting 
environments.  Radio  Sci,  vol  I 2.  1977.  p 75-87 

6 I 110  SI  I . and  Wait,  JR.  IM  Wave  propagation  in  a laterally  nonumfomi  troposphere.  Cooperative 
Institute  lor  Research  in  environmental  Sciences  (OIRPS).  University  ol  Colorado.  1 M Report  I , 

June  I 5. 1977 

7 Hamming.  KW  Numerical  Methods  for  Scientists  and  engineers.  McGraw-Hill  (Second  edition).  McGiaw 
Hill.  New  York.  1973 

X Mori  it  t.  DG.  and  Shell  man.  Oil,  "MODI  SR(  II  an  impioved  computer  program  I oi  obi  am  mg  I I I Vll  II 
mode  constants  in  an  earth-ionosphere  waveguide.  Naval  electronics  Laboiatory  Center  interim  report 
771  prepared  lor  the  Detense  Nuclear  Agency  (DNA).  I October  1976 


any  phy  sk  .ill  \ impo:  I .in  I tic  l.mgulai  region  ol  1 1 u-  complex  o i 1 1 v . 1 1 1 1 sp.k  c In  I In  s repot  1 
Ilk-  lit  ill  l > ol  llu' KIM  lot  finding  modes  associated  kith  tropospheric  duels  over  the  Ire- 
1 1 ik* in  v t.mgc  liom  \ III  to  S Kuut  kill  Ik'  examined  ( onip.it  isons.  pm ticulurly  in  rcuard  to 
compute!  linn.'  .mil  i'.im'  o|  inipli'iiu'iil.ilion.  kill  he  in.uK'  with  lln.'  more  conventional 
V'kl on -R.ipllxon  si  Ik'llk'  w Ilk'll  1 1 .is  In l ilk'll  ilk'  Ikisis  lor  Ilk’  pivv lolls  lloposplk'l  k wave- 
guide i .iK  til.it  ions  prose  ni  oil  in  rolo  roik  os  I through  (>  I ho  so  oonip.insons  kill  ho  made  lor 
•i  ground  based  iluol  produced  In  .in  olov.iloil  l.iyor  holwoon  .ihoul  iK.'in  .mil  305m  wlnoli  is 
characterized  In  ,i  relr.ictivity  deficit  ol  about  40.  I ho  I. nor  si  met  mo  appropriate  lo  this 
ground-based  duct  is  .ipproMin.itoil  In  .i  Irilmo.ir  niodol.  In  onoli  lino.ir  l.iyor  I ho  h.isk  field 
solulions  .iro  expressible  in  t onus  ol  modi  lied  ll.inkol  functions  ol  ordor  ono- third  .ind  tlioir 
derivatives  as  defined  In  the  si. ill  ol  the  ( omput.ition  Laboratory  <f  Ilicsc  liiiiotions.no 
lino. nl\  rol.itod  to  \iry  lunctions,  li  should  ho  .ipprooi.ilod  lh.it  the  HI  M developed  in 
reloreneo  N is  ro.illy  quite  a general  complex  root  finder  and  that  in  the  present  study  it  is 
being  applied  to  a very  specialized  problem  area. 

In  the  follow' mg  section  notation  lor  some  ol  the  more  important  sy  mhols  is  given 
In  section  3 some  results  ol  complex  variable  theory  . useful  for  an  understanding  ol  phase 
contours  and  their  general  behavior,  are  presented  In  section  4 an  abbreviated  description 
ol  the  Rl  M is  given  I oi  a more  complete  discussion  ol  the  latter,  the  interested  reader  is 
relerred  to  reference  X.  \ summary  ol  the  mode  equation,  along  with  modifications 
demanded  by  the  Rl  \1 . will  be  found  ill  section  5 Appendix  \ is  a mathematical  supple- 
ment to  section  5 Results  are  presented  in  section  <>  and  conclusions  briefly  summarized  m 
si' 1 1 ion 


'*  ( i imputation  I ahoralorv . I allies  ol  the  Modified  llankol  f millions  ol  Oulu  One  I liird  and  I heir 
Directives.  II.ii  vaid  In  ner  si  tv  Press.  ( amhi  id  ge . Massachusetts,  I '05 
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2 NOTATION 


I i pure  1 shows  .1  typical  elevated  layer  structure  with  altitude  as  ordinate  and 

<2n  I i as  abscissa  with  n the  modified  retractive  index  tie.  earth  curvature  included)  I he 
linear  regions  are  denoted  by  1 through  3 and  the  trihnear  model  serves  as  the  basis  for  the 
\tud\  reported  here  NIKS  units. ire  used  throughout  the  study  and  some  ol  the  more  im- 
portant symbols  used  are  the  following: 

x.  y . / Cartesian  coordinates;  / is  measured  vertically  upwards,  and  x is  meas- 
ured horizontally  in  the  plane  of  incidence 

u)  2rr  X frequency  of  rt  wave 

t Original  modal  equation 

I ■’  i = 1 . 2.  3.  4 represents  pole-tree  discontinuous  version  of  modal 

equation 

M i=l.2.  3.  4 represents  pole-free  continuous  version  of  modal  equation 

R|j  Plane  wave  reflection  coefficient  from  everything  above  height  1)  in  the 

due  t 

R|j  Plane  wave  reflection  coefficient  from  everything  below  height  I)  in  the 

duct 

0 Complement  of  the  angle  of  incidence 

s sin  0 

c cos  l ) 

k f ree  space  wave  number 

II  Height  where  modified  index  of  refraction  is  taken  equal  to  unity  II 
will  always  be  taken  in  region  I (see  figure  I ) 

Pj  iNi  Argument  of  llankel  function  at  lower  boundary  of  region  N (see 

figure  I l 

p-)(Ni  Argument  of  llankel  function  at  upper  boundary  of  region  N (see 

figure  1 l 


?<n  II 


I igure  I Schematic  ul  trihnear  model  (n  is  the  modified  refractive  index  referenced  to  unity  at  the  ground). 
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\igument  ol  I l.mkel  function  .it  lioiylil  I)  in  tegnm  1 

i 12'  Slopes  detmed  in  ligure  I such  that  o j mil  o ^ (l.iiulos 
< 0 

x^ 

1 l.mk s‘l  liiiKtioiisol  or de i one  third  .is  iletmed  In  * ompulation 
I ahoralorv  1 

Derivatives  ol  ll.mkel  luiK  lions  ol  order  one  thud  .is  delmed  In  < om 
put  .it  ion  I ,d'i iratorv  1 

I KMiel  idle.  turn  coelficient  lot  the  ground  lot  I I pol.n i/.itioti 
surta.e  i ms  hump  height 
( iround  «. undue livitv 
( in xin d permit  ti\it\ 

I ree  sp.tee  permittmt> 

1 1 me 

\ constant  equal  to  0 X 5 3(>f > ^ 2 I KS3H‘>s  I 
Imaginarv  part  ol 
Real  part  ol 

I unction  ol  a complex  variable  win. Ii  is  .male  tie  everv  where  except  . 
a finite  number  ol  isolated  poles 
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I or  understanding  much  ot  the  material  presented  in  tins  report  it  is  desirable  that 
the  readei  have  at  least  a rudimentary  understanding  ot  phase  contours  and  their  general  be- 
havior tor  functions  ol  a complex  variable  which  are  analytic  everywhere  except  at  isolated 
poles  I his  section  is  intended  to  serve  that  purpose  I he  objective  is  not  to  present  proofs 
or  theorems,  but  to  imparl  briefly  the  requisite  bac  kground  tor  understanding  the  topology 
ol  pll  ase  contours 

I el  ( a vvl  be  a function  of  a complex  variable  which  is  analy  tic  every  where  except  at 
isolated  poles  I he  phase  of  (it  w)  is  defined  as  arg  (it  vv  I tan- 1 ImKiiRet(i)  With  allow- 
ance for  the  sign  ot  Imlti  I and  RelG  1 this  phase  is  single  valued  in  the  interval  0 ' ' arg(. 

■ .'(it)  except  at  a /ero  or  pole  of  the  function  where  the  phase  is  undetermined  I xpansion 
ot  such  a tunction  about  any  simple  zero.  wq.  ol  the  function  provides  some  insight  into  the 
phase  behavior  of  the  function.  Assuming  that  w is  close  to  wy.  expansion  to  lust  order  is 
sufficient  and 

( it  VV  I ( it  W y ) + ( i(  W ) ( W - Wy  I 

( ( Wy)  e‘^0  |Aw,  e‘^  . Aw  w - Wy 

(,'twyi  Aw  e"^  + <?l  (I  . 

Polar  notation  has  been  employed  and  do  is  a constant,  namely  , the  phase  of  (it  wl  at  wy 
while  0 is  the  phase  ot  w - wy  as  illustrated  in  figure  2(a).  I his  indicates  that  near  a simple 
/ero  ol  ( a w ) the  phase  ol  (it  wl  varies  linearly  with  p and  thus  increases  in  a counterclock- 
wise sense  about  the  /ero.  I his  leads  to  the  phase  diagram  in  figure  2(b).  \ set  of  lines  ol 

constant  phase  thereafter  referred  to  as  phase  contours)  ranging  from  0 to  In  radians  ema- 
nates radially  (solid  lines)  from  a simple  zero  I lie  dashed  lines  depict  possible  phase  con- 
tour behavior,  in  the  region  beyond  the  neighborhood  of  wy.  in  order  to  emphasize  that  m 
this  region  the  phase  contours  are  generally  not  radial.  In  view  ol  the  phase  behavior  near  a 
zero  ot  (itw)  it  is  conceptually  useful  to  define  a zero  ol  (i(w)  as  a “source’' of  a set  ot 
phase  contours.  I or  an  ti*'1  order  zero  the  expansion  ol  (i(w)  leads  to  similar  phase  behav 
lor  m the  neighborhood  of  the  zero  with  the  exception  that  the  phase  of  (it  w ) varies  linear- 
ly with  up  It  is  also  true  that  the  phase  change  of  (it  w)  around  a counterclockwise  contour 
enclosing  an  n^1  order  zero  of  (it  w)  is  2nrr 

I he  phase  behavior  of  fit  w)  near  a pole,  Wp.  ol  (it  w)  can  be  obtained  by  examining 
the  Laurent  expansion  ol  (itw)  about  Wp.  Using  a simple  pole  as  an  example,  the  Laurent 
series  is  given  by 


c»o 

(itw)  = N an  (w  - wp)n  . 
n = -I 

Assuming  that  w is  close  to  Wp.  a single  term  is  sufficient  and 
a_  i i0i 

G(  w)  = w _ w = la_  i|e  I A w P ^ e-1^*  • Aw  vv  - w„ 

p 1 P 


- I 


Aw 


,itd- 1 -0i 
e 


(2) 
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I idlin'  - I’liaso  con  tout  heli.ivioi  neat  .1  /cm  ol  Uwl 

I’h.icc  com  1 on  1 in  (lie  ncigliboihooil  ol  w0 

Phase  con  1 0111  lies  oml  1 lie  iieiplihoihoiul  ol  wo 


w he  re  | is  simply  .1  constant  (a  Laurent  expansion  coefficient  l.  | is  also  a constant 
equal  to  the  phase  ol  a_| . and  0 is  the  phase  ol  w - Wp.  Unis,  near  a simple  pole  ot  (i(wl 
the  phase  of  (it  w)  varies  as  -0  <-n0  for  an  11"'  order  pole  I This  leads  to  a phase  behavior 
similar  to  that  for  a zero  of  (it  w),  except  now  the  phase  ot  (it  w ) increases  in  .1  clock  wise 
sense  about  the  pole  as  illustrated  in  figure  3.  for  reasons  which  will  become  apparent  in 
the  next  section  a pole  is  defined  to  be  a "sink"  ol  a set  ol  phase  contours. 

A summary  ol  the  basic  principles  1 list  dealt  with  is  as  follows  for  a (unction  ol  a 
complex  variable  which  is  analy  tic  every  where  except  at  isolated  poles 

I Phase  contours  do  not  intersect  except  at  zeros  or  poles 

2,  Zeros  and  poles  act  as  “sources"  and  "sinks,"  respectively  , ol  phase  contours, 
and  the  phase  behavior  is  radial  in  the  neighborhood  ot  either  a pole  or  a /ero. 

3 Near  a zero  the  phase  increases  in  a counterclockwise  sense  about  the  /ero. 

4 Near  a pole  the  phase  increases  in  a clockwise  sense  about  the  pole 

1 lie  above  principles  lead  to  the  conclusion  that  the  phase  contours  lor  a I unction 
with  a single  /ero  or  pole  tin  the  finite  region  ol  the  complex  plane  I are  simply  radials  ex- 
tending from  the  zero  or  the  pole  to  infinity.  Ol  more  practical  interest  is  the  phase  behav- 
ior associated  with  a function  ol  many  zeros,  of  many  poles,  or  ol  a mixture  ol  the  two. 

I he  phase  contour  behavior  for  such  functions  is  dictated  by  two  general  characteristics 
One  characteristic  is  the  general  phase  contour  behavior  associated  with  a neighboring 
pole  and  zero.  I he  other  characteristic  is  the  general  phase  contour  behavior  associated 
with  two  neighboring  zeros  or  two  neighboring  poles  I liese  characteristics  are  illustrated 
by  ligure  4 f igure  4(a)  shows  phase  contours  for  a function  containing  one  simple  zero 
and  one  simple  pole.  The  figure  illustrates  the  tendency  for  the  combination  ol  a zero  and  a 
pole  to  act  as  a "source"  and  a “sink"  of  phase  contours  (assignment  ol  "source"  to  a zero 
and  "sink"  to  a pole  is  arbitrary  and  the  reverse  assignment  would  be  equally  descriptive  tor 
the  purposes  intended  here)  A phase  contour  in  such  a situation  is  terminated  on  one  end 
by  a zero  and  on  the  other  end  bv  a pole.  I Ins  termination  ol  the  phase  contour  in  Ihe 
I mite  region  ol  the  complex  plane  is  referred  to  as  "phase  contour  extinction  ' I igure  4 ( h 1 
shows  phase  contours  for  a function  of  two  simple  zeros.  I he  phase  contour  behavior  is  in 
keeping  with  the  idea  that  the  two  zeros  are  both  sources  ol  phase  contours.  A similar  situa- 
tion results  lor  two  poles  with  the  exception  that  the  phase  contours  have  an  opposite  sense 
o|  rotation  about  the  poles  I lie  concept  ol  "source"-‘Nink"  pairs  is  not  violated  here  as  the 
phase  contours  are  terminated  by  the  appropriate  "sinks"  or  “sources"  at  inlimty  One 
might  wonder  il  the  phase  contour  associated  with  I wo  zeros  t or  poles  I might  ever  exhibit 
"phase  contour  extinction  ” I lie  suggested  behavior  would  be  that  shown  111  ligure  4(a  I 
with  the  pole  (or  zero)  replaced  by  a zero  (or  pole).  It  can  be  seen  why  such  a case  could 
never  occur  as  the  sense  of  rotation  ol  the  phase  contours  about  one  ol  the  zeros  (or  poles) 
would  be  contrary  to  the  precepts  put  forth  in  principles  3 and  4 above,  figure  4 illustrates 
phase  contour  behavior  for  relatively  simple  functions.  However,  as  will  be  seen  later  the 
phase  behavior  ot  the  more  complicated  functions  of  concern  in  this  study  possesses  Ihe  same 
basis  characteristics  as  those  described  above. 
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\ I ho  irem  i'I  1 1 mi  pit.'  x v ar  table  theory  states  that  lor  a I it  net  ion.  Ot  w t.  ol  a v <>m 
pU  \ \anahlc  with  isolated  poles  hut  otherwise  analy  tie.  the  accumulated  phase  , hangc  ot 
l pe  tmu  tioa.  around  a closed  counlerclov  kwtse  eontoui.  ( m the  eomples  plane  ol  the 
tmietioii  argument . is  equal  to  2n  times  the  numhei  ol  zeros  munis  the  number  "I  poles 
eneloseil  h\  ( where  n’1'  order  zeros  and  poles  are  counted  n times  dip  si  In  addition  no 
zeios  ot  poles  mat  exist  on  ( I he  interested  reader  is  reierted  to  lelereike  K)  toi  tin 
pioot  ol  this  theorem  In  passing  it  is  pointed  out  that  this  theorem  suggests  .ill  tin  mn 
. epts  disi  ussed  m the  prev ions  seetion.  I he  theorem  nulls ates  that  the  phase  , hange  ol  .1 
turn  1 1 mi.  < ■(  w l.  .ihout  a elosed  counters  hu  kw  ise  eontoui .(.  cm  losing  only  a set  ol  n sim 
pie  zeios.  is  equal  to  2ii7T  I Ills  directly  implies  that  every  phase  eontoui  assos  lated  w itli 
e.u  h zeio  crosses  ( I lie  Kl  M was  developed  Irani  this  basis  pnmiple 

Some  fundamentals  ot  the  ap  plica  turn  ot  the  Kl  M lot  ind  zeros  ol  a I unci  ion.  ( it  w t. 
illustrated  in  figure  <>  \ search  res  tangle  is  placed  about  some  region  ol  the  complex 

pi. ms  I he  ss'.u s h res  tangle  is  sli\  icleil  into  a grid  ol  mesh  squares  whose  corners  w ill  be 
,..,||ed  mesh  points.  I he  mesh  square  size  is  optional  and  is  usually  selected  according  to  the 
expes  ted  zero  spacing.  \ mesh  size  that  is  too  huge  can  result  in  complications  which  max 
nltnnatelv  lead  to  missed  zeros  as  will  be  explained  m section  (>.  However,  a mesh  size  that 
|S  |,,n  small  leads  to  a costly  extended  computer  run  time.  Keginning  at  the  upper  lei t coi 
net  ol  the  sears  It  re.  tangle,  a boundary  search  tor  I)  or  I SO  phase  contours  is  conducted  m 
the  counterclockwise  direction  \ny  phase  eontoui  would  do:  however,  the  0 and  ISO 
phase  contours  are  selected  because  mathematically  they  are  easily  located,  occurring  when 
| in  l ( i ) 0 I he  search  is  conducted  by  evaluating  (itw)  at  the  mesh  points  along  the  search 

IC(  I. ingle  boundary  \V  lieu  lint  < > I s lianges  sign,  it  indicates  that  a 0 or  I SO  phase  contour 
has  just  been  passed  (points  V I).  and  (,i  Once  either  ot  these  phase  contours  is  located, 
the  bouiulai y search  is  temporarily  halted  while  the  0 or  ISO  phase  eontoui  is  traced  into 
the  mtei  nu  ol  the  seal  eh  res  tangle  by  inspection  ot  ImtOl  at  the  corners  ol  the  mesh 
squ.iics  leountcu  lockwise  inspection  beginning  at  the  top  let l corner  ol  each  mesh  square! 

| he  phase  i ontoui  is  followed  either  until  a zero  ol  ( it  wl  is  discovered  ( points  K and  I I or 
until  the  seal,  h re  t angle  bound  a i y is  encountered  (as  would  be  the  case  lot  the  phase  e on 
lour  between  ( . and  III.  one  ot  which  will  always  occur  provided  no  poles  exist  in  the  interi- 
or ot  the  search  res  tangle  \\  hen  a zero  is  located  (by  means  which  will  be  explained  short 
ly  I.  its  losatiou  is  saved  I hen  the  phase  contour  is  traced  out  the  opposite  side  ol  the  zero. 

Ill  (Inn,  lull  K\  Kiwi  IW  and  Vim  lev  Kl  < , implex  V.n  i.iblcs  and  \pphc.il  inns  | I Inrd  edilinnl. 

\te(. raw  Hill  Vw  Vmk.  1974.  p M(, 
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ig mo  (>  Root  finding  method  ( Rl  M I lor  a function  (r<w). 

Phase  contours  for  (>(wl 

V)  Zeros  o|  (,<  wi 

having  undergone  a I HO  phase  change  (see  tig  2(b)),  until  the  search  rectangle  boundary  is 
again  encountered  (points  ( and  I I.  When  the  phase  contour  exits  the  search  boundary  . such 
as  at  points  ('.  F.  or  II.  the  mesh  square  which  contains  this  occurrence  is  (lagged  so  as  to 
avoid  following  that  particular  phase  contour  again  at  a later  time  during  the  boundary  search 
Also  at  such  a point  (point  C.  F or  II)  the  phase  contour  trace  is  stopped  and  the  boundary 
search  is  resumed  at  the  point  where  the  lust  0 ' or  180'  phase  line  was  encountered  (eg. 
points  A.  1).  or  G).  Once  the  entire  search  rectangle  boundary  has  been  inspected,  all  the 
zeros  ot  the  function.  G(  w).  located  within  the  search  rectangle  will  have  been  found. 

I he  location  of  a zero  is  evident  by  the  intersection  ot  phase  contours  (section  3. 
tig  2(bi)  Iheretore.  the  intersection  ot  the  0°  or  1X0  phase  contour  with  any  other  phase 
contour  locates  a zero  of  G(  w)  File  other  phase  contour  chosen  lor  t Ins  purpose  is  the  00 
or  270  phase  contour,  again  chosen  for  simplicity,  as  these  contours  are  easily  recognized, 
occurring  when  RelG)  = 0 While  a 0 or  1X0'  phase  contour  is  being  traced.  RetGl  is 
also  examined  at  the  corners  of  each  mesh  square  to  locate  a RelG)  sign  change  which  indi- 
cates that  a 00  or  270  phase  contour  has  entered  the  mesh  square  (the  bold  lined  mesh 
squares  in  tig  7(a)  and  (h(l  Such  an  occurrence  indicates  that  a zero  is  probably  within  that 
mesh  square  dig  7 (a ))  or  perhaps  within  an  adjacent  mesh  square  ( fig  7(b))  Once  a mesh 
square  is  known  to  contain  a zero,  a more  precise  location  of  the  zero  is  obtained  by  an 
interpolation  scheme  which  employs  both  the  magnitude  and  phase  ol  the  function  G(w). 
i For  details  ol  this  interpolation  scheme  the  reader  is  referred  to  reference  8.1  F ollowing 
this  a Newton-Ruphson  iteration  pinpoints  the  location  of  the  zero 
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I lie  KIM  torbids  the  existence  ol  poles  in  the  enclosed  region  01  on  the  boundary  ot 
the  search  rectangle.  Such  occurrences  could  lead  to  undetected  /eros  las  will  he  explained 
in  section  (')  due  to  "phase  contour  extinction"  (see  section  3.  tig  4la)l  Poles  may  exist 
external  to  the  search  rectangle,  hut  even  under  such  restrictions  it  the  poles  are  in  proximity 
to  the  search  rectangle,  they  can  be  a definite  source  ot  problems  due  to  their  eltect  on  the 
phase  contours  associated  with  zeros  that  are  interior  to  the  search  rectangle  I Ins  causes 
problems  similar  to  those  caused  by  the  increased  density  ol  phase  contours  shown  in  figure 
I 2 and  discussed  in  section  (>  Therefore,  it  is  best  to  have  a pole-free  function  in  the  finite 
part  of  the  complex  plane  in  order  to  assure  optimum  performance  ol  the  KI  M.  In  effect, 
all  poles  should  be  transferred  to  the  point  infinity  in  the  complex  plane 


5 SI  MM  AKA  (>l  V\  \\  I (,(  11)1  I (.)(  A I IONS 

In  this  section  Ilk'  nrijiin.il  mod.il  >.-< | u.i I ion  lot  tin.*  Inline, n model  ol  1 1 jin n."  I along 
with  niiulilk  .it  ions  ilk  1. 1 led  In  the  K I M will  In  iimm.in/cd 

I Ik  lund.imcnt.il  p inlik'ii i in  pi. i ii. ii  w.nce  tilde  m.iksis  is  i Ik-  solution  ol  tin.'  modal 

Cipi.ltlOII 

t (it  i I - K]  jit)  i K|j(tt  i U t s i 


for  eigenangles  Five  formulas  summarized  below  lot  the  plane  wave  reflection  cocffi 
. i.  it-.  K|>  .ms]  l<|)  .us  tot  hori/onl.il  pol.m/.ilion  ,nul  .issiime  .i  tune  ilepeiuienee  ol  c,t*;l 
I oi  i derivation  ol  the  formulas  the  interested  readei  is  referred  to  refer*  nee  I 1 [’he  re 
tl.  tion  . ost  Ik  lent  K | ) is  i .ik  ill.  ited  In  e.ile ul.it  mg  in  seiiueik  e I lie  follow  ini’  i|  n.  in  lilies 
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\lso.  the  quantities  N / 2 \A2.  NZ  I . etc.  are  the  numerators  ol  . A 2.  • | etc.  and  similarly 
l)/2.  1 ) A 2 . I)/ 1 . ete.  are  the  denominators  of  . A2.  | etc 

Beeause  the  RIM  requires  no  poles  within  the  search  rectangle  and  preferably  none 
in  the  Unite  part  ot  the  complex  plane  (see  sections  4 and  <> I.  the  modal  equation  requires 
mollification  as  discussed  subsequently  I he  significant  point  at  this  juncture  is  that  the 
modification  can  lead  to  overflow  problems  it  the  wavellelds  become  very  evanescent  or 
exponentially  small  m regions  2 or  4 (see  tig  I ) I he  overflow  problems  result  from  the 
quantities  . s and  At  2 ) defined  in  equations  (4  > and  ( 5 ).  respectively  In  such  a circum- 
stance. the  overflows  can  be  avoided  by  calculating  .|  as  follows 


Its  ( pj  ( 2 It  - c^ni  It  | ( p|  1 2 l ) 
mas,  ki1  3 (lv’o  (P|(2»-e4ffi  3 h,  (p|(2n) 


\/.l 
I )/ 1 


( I 4 l 


( fbserve  that  • s and  At  2 1 as  defined  by  equations  (4  ) and  1 5 I now  do  not  appear  explicitly 
in  equation  ( 14  i I bus  it  the  rf  wave  is  sufficiently  evanescent  in  regions  2 and  3.  equation 
i 14  i replaces  equation  HO  and  there  is  no  need  to  calculate  . s and  At  2 > via  equations  (4t 
and  ( 5 ) I his  is  an  example  of  what  is  referred  to  as  function  switching  at  various  stages  of 
this  report  As  will  be  discussed  later,  the  decision  to  use  equation  ( 14  I is  based  on  the 
magnitude  of  the  quantity  I | . which  is  defined  as  follows: 

I I = ! Re  (i  (p  | ( .’mi  ’ -) 1 + yj  | Re  (i  ( p o ( 2 ))■'  -)|  ( I 5 1 


rile  reflection  coefficient  R|)  is  calculated  by  calculating  in  sequence  the  following 
quantities 
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I he  quantity  Pj  ( I I is 

1/  s 

P|  1 I I = (k  U|  I ' (S“  -Q|  II)  (|4) 

As  before.  N/l).  NB1 . etc . denote  the  numerators  ol  .().  B(  I I.  etc,  whereas  DZO.  DBI , etc. 
denote  their  denominators 


\jt.mi,  because  ol  modifications  required  In  t lie  Kl  \l  overflow  problems  can  o>  <-  i i i 
il  1 1 10  it  fields  become  very  evanescent  « u exponentially  sm.ill  .it  the  ground  In  this  case 
mound  quantities  i .in  he  ohm  nut  eil  In  using  the  I nl  lowing  .iltorn.it  ivo  expression  lot  R j j 


M ) 


S ( ll  s I pjjl  1 II 
S (ils  (p|jl  I II 


It:  i 


' h j (p|)l  I ll)  I I ( d | k I1  ■'  (Il  s (P|,l  I II  - e l7ri  ' li’|  (P|)l  I ll) 


.■kill  3 


' ll  | <|'|)l  I ll)  ~ I tO  | k I ' (ll's  (P|)l  I 


. 1 77 1 3 


\ K 1 1 
I ) K i < 


' li’|  <l’|)l  1 ") 
(20l 


I lie  decision.  .i'  w ill  he  dis.  usseil  presently  . to  use  ei|ii.ition  ( 2()i  nisle.ul  «>l  equation  i I hi  is 
based  upon  the  magnitude  ol  the  quantity  . I s.  defined  .is  lollows 


I s -=  Re  ( i < P i i 1 1 1 ' ") 
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\ rutlier  considerable  overestini.ite  ot  the  magnitude  ol  the  modal  equation.  niodi- 
t ied  so  .is  to  remove  .ill  poles  to  the  point  infinity  in  the  complex  plane  is  given  In 


exp  1 1 


I si  . 


where  I ; is  defined  by  equation  ( 1 5 I and  Lm  In  equation  ( 2 I I I lius.  no  lunelion  switching 
is  reqiureil  it 
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whole  expi  l.i  is  the  maximum  numher  the  computer  can  handle  (exp(  I I will  he  referred  to 
as  the  il \ namie  range  ot  the  machine)  It  I | + I s I the  simplest  decision  lor  lunction 
switching  is  to  use  equation  I I4i  instead  ol  equation  H>l  it  I | > I 2 and  to  use  equation 
I 20i  instead  ol  equation  t IX)  it  J '-12  I or  reasons  of  implementation  which  will  he 
• mie  apparent  later,  it  is  a dev nled  adv .ullage  when  using  the  R I M to  avoid  t unction  sw itch- 
ing as  mm  h as  possible  In  this  sense  a sw  ik  lung  scheme  whu  h hotter  utilizes  the  dynamic 
tango  ot  the  machine  is  indicated  hy  the  tollowmg  llowchart 


Iii  the  present  study  the  Univac  1110  has  been  used  with  single  precision  format  lor 
the  Kl  M so  that  the  dynamic  range  is  on  the  order  ol  1(P*  or  ~ e*7 

\s  has  been  mentioned  previously  (also  see  sections  4 and  6).  the  Kl  M requires  that 
there  be  no  poles  within  the  boundaries  of  the  rectangular  region  being  searched  for  solu- 
tions. and  in  tact  it  is  preferable  that  all  poles  he  transferred  to  the  point  infinity  in  the 
complex  plane  Since  the  denominators  ol  the  preceding  expressions  tie,  DZ:.  DA:.  DZI 
etc  i possess  zeros,  the  elimination  of  all  poles  within  the  finite  part  of  the  complex  plane 
can  be  assured  by  multiplying  the  t function  (equation  (3))  through  by  the  denominators 
Because  ol  the  switch  conditions  based  on  the  parameter  I (see  flowchart  on  page  lb)  there 
are  iii  effect  four  distinct  representations  ot  the  modal  equation.  I he  pole-free  versions  ol 
these  are  presented  below 


Pole- free  Discontinuous  Version  of  the  Modal  liquation 

1)  1 | + I s < 1. 

i | = dz:  x da:  x dzi  x dai  x drd  x dzo  x dbi  x dkd  x t C4i 


w here  DZI  is  given  by  equation  ( (> ) and  DKD  is  given  by  equation  (IK) 
n)  L|  + Is  > I.  and  I | > l.i  and  I s < 1 

Is  D/I  X DAI  X DKD  X DZO  X DBI  X DRD  X 't 

w here  DZ  I is  given  by  equation  ( 14)  and  DKD  by  equation  ( lb) 

in)  L. | + I s > L and  I. s > L|  and  I | < l 


1 , = DZ2  X DA2  X DZI  X DAI  X DKD  X DKD  X t 

where  DZI  is  given  by  equation  (bland  DRD  by  equation  ( :() ) 
iv)  MiN  ( L, . Ls)  > L 
I 4 = DZI  X DAI  X DKD  X DKD  X ‘f 


( :.i  i 


( :o  i 


<:: 


where  DZI  is  given  by  equation  (14)  and  DKD  by  equation  ( :0 ). 

liquations  ( :4  I through  ( 27 ) are  pole-tree  representations  of  the  modal  equation  in 
different  regions  of  the  complex  plane.  It  is  clear  that  the  functions  I | . I s I },  and  I 4 arc 
not  continuous  so  that  they  are  referred  to  as  the  pole-free  discontinuous  version  ol  the 
modal  equation  It  is  not  surprising  that  a discontinuous  set  ot  functions  used  to  represent 
the  modal  equation  in  different  regions  of  the  complex  plane  can  lead  to  difficulties  when 
used  in  conjunction  with  the  KIM  which  requires  phase  continuity  111  order  to  trace  lines  ol 
constant  phase.  By  examining  the  asymptotic  behavior  ot  the  modified  llankel  functions  of 
order  one-tlnrd.  it  can  be  shown  (see  appendix)  that  within  the  significant  region  ol  the 
complex  plane  a continuous  representation  of  the  pole-free  modal  equation  is  as  given 
below 
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I qiulions  ( 2H>  through  i'll  describe  .1  pole- tree  continuous  version  ol  the  modal  ei|ua- 
tion  vv  his  li  will  be  discussed  Ini  tiler  in  sec  turn  (1  I lie  continuity  las  tors  introilueesl  1 11 
s-i | ua I ions  1 2b  i through  i'll  iiilroduee  braneli  points  alone  the  real  axis  I he  branch  points 
m I lie  sle  no  in  111a  tors  ol  eipiat  ions  I 21 ) 1 anil  I ' 1 I are  sinks  ol  phase  eon  tours ; liowev  er.  no 
coniphs  ations  oceui  as  these  equations  are  employs’ll  in  regions  well  removed  Iron)  the 
branch  points.  I he  branch  points  in  the  exponential  terms  ol  equations  1 2X)  through  i.'Ul 
san  cause  problems  due  to  phase  discontinuities  across  then  associated  brunch  cuts,  liowev 
ci.  it  has  been  I omul  that  these  discontinuities  are  not  large  enough  lin  the  region  enclosed 
In  reasonable  search  rectangles)  to  interfere  with  the  phase  contour  tracing  ol  the  Kl  M il 
the  branch  cuts  are  m the  lower  half  ol  the  complex  0 plane  and  sufficiently  divergent  from 
the  real  axis 
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<»  RESULTS 

A previous''  ease  studs  ('I  a ground-based  duet  produced  by  an  elevated  layer  (see 
tig  Si  between  IX3m  i <>()()  Iti  and  305m  ( 1000  ttl  eliaraeteri/ed  by  about  a 40m  unit  delieit 
1 1 >|  1 M i ix  reexamined  in  this  section  by  use  ol  the  HI  M As  shown  in  figure  S the  layer 
height  is  /|  and  the  lay er  structure  is  approximated  by  a trihnear  model  The  abscissa  is  the 
modified  re  tract  iv  its . M ( \l  - ( n - I I X It/’  where  n is  the  modified  refractive  index  I.  and 
the  ordinate  is  altitude  m metres  I lie  frequencies  ('5  Mil/.  I 70  Mil/,  520  Mil/,  and  3.3 
(ill/  used  in  the  previous  study  are  reexamined  so  that  a direct  comparison  ol  computer  run 
tunes  sail  be  made.  I he  mode  attenuation  cutolt  was  0.5  dB  km  lor  the  three  lowest  fre- 
quencies and  I cl B km  for  3 3 (ill/.  \s  a result  the  number  ol  modes  used  ranged  Irom  2 at 
<i5  Mil/  to  X I at  3.3  ( >11/ 


figure  8 I nlmeai  modified  reluctivity  profile  referenced  to  zero  at  the  ground 
BOUNDARY  SEARCH 


figure  l>  Application  ol  the  Kl\l  to  I lie  l'2U  Mil/  case 

0 and  180  phase  contours 

@/  Zeros 


Direction  ol  search  path 
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B\  use  ot  ilk  Kl  \|  mil  ilk'  poll  l ii\-  m ml, 1 1 Imi.  I k 111  I | ill- 1 mi  ll  by  1 1 1 i.ili on  i .'  4 I ol 

sc,  Ill'll  ' I III-  ClgCIlallglcs  lol  I IlC  llirfC  lowest  lleqUCIk  ICS  WCIC  obtained  WllllOUl  lOllipIli  ,1 

I ion.  il  ik'  iti  Ilk-  lack  oi  tune  lion  switching  < dealt  with  m \ik  lion  5 i .11  these  relativ  cly  low 

I I i'i|lli'IU'li‘s  Ik’  onls  tlk'  llllk  Ill'll  I | was  IISCill  I Iglllc '>  IS  ,i  k presentation  ol  lilt'  Kl  M 
applied  to  Ilk-  •'2D  Mil/  case,  ilns  heme  an  interesting  represeiilatise  ol  1 1 1 is  lower  frequency 
domain  I lie  ordinate  is  Iiiiid  I and  the  absc  ixs.i  is  KciW  i I lie  axes  are  scaled  dillerently  in 
oidei  lo  piesent  a clear  and  complete  siess  ol  the  essential  phase  contour  behavior  ol  the 
modal  link  lion.  I | . m the  cut  ire  region  ol  the  seau  li  rec  tangle  I sselve  /eios  lot  modal 
solutions)  are  indicated  m the  region  ot  theta  space  shown  I liese  modal  solutions  are 
presented  m table  I Since  a general  desciiption  ol  the  Kl  M lias  been  given  previously 
is.-ctiondi  only  certain  interesting  aspects  ot  figure  will  be  singled  out  I lie  rectangle 
bou ndar \ search  begins  at  the  uppei  lei  I corner  ot  I he  scare  h rectangle  and  pioceeds 

as  indicated  by  the  airovvs  In  the  diagram  \ depicts  an  example  ol  a point  vvlicie  the 
li.icuiL’  ot  ,i  phase  contour  stops,  having  encountered  the  search  rectangle  boundary  lliis 
location  is  flagged  so  as  to  av  >ul  the  i c- 1 raving  ol  that  phase  conloui  at  a later  tune  during 
lb.  boundary  search.  Il  is  interesting  to  note  that  I he  inodes  are  not  found  m u msec  ulive 
oidei  Due  lo  the  phase  contour  configuration,  the  ninth  mode  is  louiul  helore  the  eighth, 
which  m turn  is  not  lov  a ted  until  last  As  i an  be  seen  Irom  the  I igure.  the  uc  tual  < irdet  in 
which  the  modes  arc  bniikl  is  I through''  12  II  It).  S I w;elve  modes,  consistent  with 
the  previous  work  irvl  5 l have  been  lound:  however,  unlike  the  Inklings  ot  the  previous 
vv  oi  k w hk  h i elied  on  the  New  ton-K.iphson  iteration  it  is  now  cei  lam  that  all  modes 
within  the  scare  h rectangle  have  been  round 

Phase  contours  for  the  original  modal  ei|U.ilion  toll  I - RK  (ei|ualion  ' i are 
shown  in  l igure  Id  lor  a small  region  ol  the  .oniplex  plane  lor  the  !rei|ucncv  a a (ill/  I he 
abscissa  is  Kell)  I while  the  ordinate  is  lino))  I he  modal  equation  contains  poles,  lour  ol 
which  aic  indicated  in  the  ligure  along  with  six  zeros  I he  zero  at  the  origin  is  a mode 
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1 igure  10  Phase  contours  tor  the  iMigm.il  modal  equation 
(9>  Zeros 
* Poles 


winch  has  zero  excitation,  therefore  it  is  o|  no  concern  I lie  ellect  ol  the  poles  on  the 
phase  contours  is  clearly  evident  in  the  complicated  phase  contour  structure  “Phase  line 
extinction'  and  its  disastrous  results  can  he  seen  lor  the  hy  pothetical  search  rectangle 
shown  (dashed  line!  I lie  problem  is  that  not  all  the  phase  contours,  associated  with  the 
zero  at  Re(tf  I s 2 (>  X 10“-’  rad.  cross  the  search  rectangle  boundary  some  are  confined  to 
a limited  region  in  the  neighborhood  ol  the  paired  pole  and  zero.  As  can  he  seen,  the  0 and 
I K(J  phase  contours  ol  that  zero  never  cross  the  search  rectangle  boundary  therefore  the 
zero  would  not  he  located  by  the  RIM  This.  then,  clearly  illustrates  the  necessity  o!  a 
pole-tree  modal  equation,  at  least  within  the  interior  of  the  search  rectangle,  to  assure  no 
missed  modal  solutions. 

I igure  I I shows  phase  contours,  in  the  same  region  ol  the  complex  plane  as  t igure 
10.  lor  the  modified  modal  equation  which  is  pole-tree  and  continuous  as  discussed  m sec- 
tion 5 ( see  in  particular  equation  ( 3 I I and  the  associated  discussion  i In  contrast  to  the  case 
lor  the  original  modal  equation  (fig  Mil  the  absence  ot  poles  results  m a much  more  relaxed 
and  simplified  phase  contour  behavior  I or  such  a situation  the  ()  and  ) SO  phase  contours 
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figure  I I Phase  contours  lor  the  pole-tree  continuous  version  ol  the  modal  equation  at  3 3 tillz 
(it  Zeros 
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I igure  13  Phase  contours  in  the  region  ol  a highet-order  mode  lot  two  versions 
of  the  modal  equation  at  3 3 GH/ 

• Zero 

(7)  Pole-free  continuous  version  ol  the  modal  equation 
© Pole-tree  discontinuous  version  of  the  modal  equation 


responsible  lor  the  spurious  phase  contour  behavior  is  illustrated  In  figure  13  f igure  13(a) 
shows  phase  contours,  in  the  region  of  a single  mode,  lor  the  pole-tree  continuous  modal 
equation  (equation  ( 2Ki  of  section  3 ).  f igure  13(b)  shows  phase  contours  in  the  same 
region  for  the  pole-free  discontinuous  modal  equation  (equation  ( 24  i of  section  5 ) As  can 
be  seen,  the  absence  ol  the  continuity  factors  in  the  discontinuous  vetsion  removes  the 
presence  of  the  spurious  phase  lines  I he  resulting  simplified  phase  contour  structure  allow 
a significant  reduction  in  computer  run  time. 

I he  temptation  to  simply  use  the  discontinuous  version  must  be  approached  with 
aii lion  I igure  I 3 illustrates  a danger  of  using  discontinuous  f unction  with  the  KIM  I he 
same  regions  ol  the  complex  plane  are  shown  in  (a)  and  (b)  of  the  figure  In  particular, 
phase  contours  ol  I | (equation  ( 24)  of  section  3 ) and  I s (equation  ( 23  ) of  section  3 I are 
plotted  in  figure  13(a)  and  figure  13(h).  respectively.  I.aeh  function  is  used  in  a separate 
region  of  the  complex  plane.  I he  dashed  line  denotes  a switch  boundary  I-  j is  used  in  the 
region  to  the  right  ol  the  switch  boundary  and  I 2 in  the  region  to  the  lef  t.  I he  type  of 
situation  presented  in  the  figure  leads  to  an  undetected  mode  as  follows.  In  the  lower  left 
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I igure  I S A problem  ol  I unction  switching  lor  the  pole  tree  discontinuous  version  ol  the  modal  equation 
Zeros 

♦ Direction  ot  search  path 
Switch  houndat) 


cornet  ot  tigurc  I 5 ( b I it  in  seen  ih.it  i s would  lx-  used  m tins  region  during  I lie  roe  tangle 
hound. u\  scutch  portion  ol  t ho  RIM  I ho  lust  <>l  tlio  lltroo  modes  in  llio  diagram  would  ho 
touiul  \ i.i  i ho  I s ]S()  ph.iso  oontoui  associated  vn  till  it  Whon  tho  hound.  iry  search  reached 
\ III.  nw  IK  h to  I | would  ho  ni.ulo,  ,ind  it  o.m  ho  soon  ill.  it  tho  (I  pluiso  c<  nit  our  ol  I > •mil 
tho  I M)  ph.iso  oontoui  ol  I | w Inch  arc  associated  with  the  second  mode  would  not  he 
onoounloiod  .ilong  tho  sonic  h p.illi  < ontmumg  on.  tho  thud  modo  would  ho  loo.itod  v la  tho 

0 pit  .in'  i on  I on  l ol  I ! \\  lion  i.  01111110  hue  k .dong  t Ik'  lop  oil  go  ot  tho  no. 11  oh  roe  tangle  hound 

.11  \ . ilk  nw  ik  h h.u  k to  I n w ou lil  ho  n i.id 0 .it  H It  in  nooii  tlr.it  tho  ()  pli.ioo  oontoui  ol  I | .1  ml 
tho  I nii  ph.iNo  oontoui  ot  I ' w Inch  lend  to  tho  noooihI  11101I0  would  not  ho  encountered 

1 horotoio.  tho  Noooiid  modo  in  tho  d1.1gr.1111  would  go  undetected  I h in  shows.  r.ithor  clear- 
ly .1  c tueial  dll  tic  tilt y oik  ( hi  11  to  rod  in  indiscriminately  employ  mg  tho  diNoontinuouN  version 
ot  tho  inoil.il  oi|u.ition  with  tho  Rl  M 

I horo  in  .1  practical  way  to  uno  tho  polo  free  discontinuous  vorsion  ol  tho  modal 
o. pi.it ion  1 ho  origin. d scale  h roe t.niglo  is  div ided  into  two  Nin.dlor  no.i roll  roo  tangles  .it  tho 
switch  boundary  Now  since  .1  single  function  is  used  in  each  rectangular  region,  tho  applica- 
tion ol  the  Rl  M to  each  separate  rectangle  will  definitely  locate  all  tho  modes. 

V tiiully  lot  tho  v ' ( .11/  case  there  wore  two  switch  boundaries  involved  ( tho  other 
hound. ns  Occurrod.it  Retd  I d-1  X lir'i.nl)  However,  instead  ot  div  tiling  the  original 
sc.  1 1 . 1 1 10.  tangle  into  throe  sinal let  scare h root. ingles  it  was  1 0111  ill  to  he  desirable  to  dis ido 
tho  original  scare  h rci  tangle  mils  at  the  second  switch  ho  tin  dare  ( Ro  1 0 ) " (>  2 X It)  ' rad  1 
and  to  use  tho  continuous  link  turns  in  tho  region  containing  tho  lust  swile  It  boundary  1 at 
Roi  m 4 "T  / I (r-'  rad  1 Alt  hough  tho  continuous  functions.  M4  and  \h  (oipi.it  ions  ( 4 I I 
and  t2'Hi  are  used  in  tho  first  search  rectangle,  their  associated  spurious  phase  lino  behavior 
is  not  seriously  dotrimont.il  to  tho  computer  run  time  (actually  M4  contains  no  spurious 
phase  behavior  while  \| ' has  such  behavior  hut  the  behavior  is  loss  prominent  at  tho  lower 
Ro  1 II 1 values ) I Ills  was  tho  procod  lire  act  ually  employ  oil  m using  tho  discontinuous  version 
ot  tho  modal  captation  with  the  Rl  M 

In  general  tor  inline. u profiles  ol  tho  typo  depicted  in  ligure  N there  would  bo  tom 
possible  switch  boundaries  and  1 01  mult  isogmon  toil  pi  ol  lies  ( more  than  throe  segments  l \ 01  y 
likolv  considerable  itioic  than  lour  It  is  cloai  that  locating  those  hound  at tes  1 01  tho  put  pose 
ot  search  1 oc  tangle  const  1 not  1011  can  present  a lot  mid. 1 bio  problem  I hero  I or  0.  tho  use  ol 
discontinuous  lunctioiis  although  possible  m principle'  can  lor  general  cases  be  very  dillicuit 
to  implement 

I able  2 lists  the  compute!  tun  tunes  lor  wav oguulc  tuns  lot  tho  two  v visions  ol  the 
R I M along  with  tho  pi ov  1011  sly  used  Newton- R a plison  method  I 01  tho  ,v -( >1 1/  case  tho 
disastioiis  ol  loot  ol  the  spurious  phase  contours  in  the  continuous  vorsion  ol  the  modal 
epilation  is  obvious  from  the  20-minute  run  time  which  is  ovei  throe  times  greater  than  that 
1 01  the  discontinuous  v orsion  ol  tho  modal  equal  ion  I ho  Now  ton  Ra plison  case  was  run  on 
an  I UM  4<>0  computer  while  the  Rl  M cases  wore  run  on  a t nivac  lilt)  computer  A rough 
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comparison  between  the  two  computers  is  that  the  I'nivac  1 1 10  is  twice  as  last  as  the 
|H\t  .too  Iheietore.  the  New ton-Kuphson  run  times  presentee)  in  the  table  are  approximate 
ami  were  obtained  simple  by  dividing  the  IHM  3t»0  run  time  In  two.  1 he  520  Mil/  case  has 
been  listed  as  an  example  ol  the  lower  frequencies.  In  this  domain  (as  a result  ol  a decrease 
m the  d\  nanuc  range  requirements!  function  switching  does  not  occur  and  spurious  phase 
line  behavior  is  not  as  prominent.  Iheietore.  only  one  rectangle  is  needed  toi  the  KIM  dis- 
continuous version  and  extremely  small  mesh  squares  are  not  necessary  lor  the  Kl  M con- 
tinuous version  1 lie  result  is  equivalent  run  times  for  the  two  methods. 

It  is  clear  from  the  preceding  discussion  that  the  KIM  is  most  easily  implemented  it 
function  switching  can  altogether  be  avoided,  lot  this  reason  a set  o!  tables  (tables  3(  a I 
through  3<  p 1 1 has  been  prepared  as  a guide  to  the  waveguide  user.  I he  tables  show  an 
estimated  maximum  frequency.  f MAX.  below  which  function  switching  should  not  occur 
toi  several  tnlinear  laser  models  and  various  search  rectangles.  I he  method  ol  calculation 
(based  on  equation  23  ol  section  5)  should  ensure  a rather  conservative  estimate  ot  I- MAX 
t te.  switching  should  not  occur  until  Irequencies  somewhat  higher  than  1MAX  are 
encountered  >. 

| MAX  is  dependent  upon  profile  parameters,  search  rectangle  location,  and  ol  course 
the  maximum  magnitude  allowed  by  the  computer  used.  All  calculations  are  lor  a machine 
, apable  * >1  handling  a maximum  mini  bet  of  e (ie,  L = 700  according  to  the  notation  ot 
section  5)  Values  ot  L other  than  700  will  be  discussed  shortly  In  all  instances  the  layer 
gradient  Q|  (see  tig  I » lias  been  taken  as  the  value  for  the  standard  atmosphere.  2.36  X 
i I0_7  nr1  Sixteen  tables  have  been  generated,  each  table  being  associated  with  a particular 

combination  ot  lour  different  layer  heights.  /[  (see  tig  Hi  and  tour  different  search  rectangles 
( the  layer  height,  / \ . and  top  right  corner.  Ret  0 l and  lint 0 ).  ol  the  search  rectangle  are  listed 
at  the  head  ol  each  tablet,  four  different  search  rectangles  are  used,  f.ach  search  rectangle 
has  a maximum  Re(0  t value  ol  1 .5  X lu_-  radians  while  the  maximum  I nit  (3 1 takes  on  tour 
different  values.  I ach  table  includes  protiles  which  have  various  layer  strengths  or  M deficits 
(1)1  LM  l ( see  tig  X i and  v arious  lay  er  gradients  a 2 and  a}  ( see  tig  I > In  addition  to  1M AX 
the  tables  give  the  maximum  possible  attenuation.  A I Ml'.  associated  with  the  search 
rectangle  ( te.  the  attenuation  associated  with  the  upper  right  corner  ol  the  search  rectangle) 

It  was  stated  that  the  tables  were  prepared  lor  a computer  that  handles  a maximum 
number  ot  e'-  where  L 700.  f or  a computer  allowing  a maximum  value  ol  L.  other  than 
700.  the  conversion  ol  f MAX  m the  tables  to  the  proper  value  is  a simple  task  as  L and 
I M AX  are  directly  proportional.  In  other  words,  it  e f is  the  maximum  value  allowed  by  the 
computer,  then  the  correct  maximum  frequency.  I(.  is  obtained  by 

t,  =—  I MAX.  <321 

1 700 

where  I MAX  is  the  tabulated  frequency 

As  an  example  ol  the  use  ot  the  tables,  we  estimate  f MAX  toi  the  present  case  study 
I tie  prot ile  parameters  and  search  rectangle  m the  table  which  are  c loses!  to  those  used  m the 
actual  study  are  /|  200  in.  Dl  LM  fO.os  -.(>  X I0_<1  in-1 . orj  .(>  X 1 0_(' m" 1 . 

Kelt!)  15/  Hr' . and  Imi  0 1 I5X  Itr'iadian  I liese  entries  are  found  in  table  3t  a), 
where  I MAX  is  seen  to  be  I 331  X It)"’  11/  I he  Univac  single  precision  value  tor  L is  87.5 
lie  let  l 87  5 in  equation  ( 32)).  therefore  no  (unction  switching  occurs  in  the  present  case 
lor  frequencies  less  than  about  ( 87  5 7(K)|  X I 331  / l()'(lH/  I 7X  It)1  11/ 
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I Ik-  tabulated  values  of  I MA  X haw  been  obtained  In  usmi: 


I- MAX 


I 


i.  I'*1"  I 1 

Kol-J - Or . / I 

\.V  a | 11 


♦ — )( s-  - 2 X l(r(’  1)1  I Al) 


3 : 


(33 1 


It  should  he  noted  th.it  the  value  ol  I MAX  is  detennmed  In  computing  the  above  quantiH 
(equation  ( 33 )l  at  both  the  upper  let t ami  upper  right  eorners  ol  the  search  reel. male  and 
selet  In il’  the  smaller  <>l  the  two  values 


I AH  I I .?<;n 


1 = 700. 

q , = 23ft  X I0"(' m'1 

Kc  1 1*  t 

15  X 10“ 1 K.ni 

InHtft  = 15  X 10"'  R;ul 

/'.  | = 200  in 

DC  - 

,h/. 

AIN 

IMA*  H /) 

A T NO 

1BA--  i 

IP  . 0 

-. 100-008 

. 1PP-P06 

1 .581*010 

19.74 

20.0 

-. 100-005 

. 100-006 

1.204*010 

15.03 

IP  . 0 

100-006 

.236-PP6 

1 . 889*0 1 0 

23.50 

29 . 0 

-. 100-005 

. 2 16-006 

1.853*010 

23.14 

IP.  P 

- . 1 P0-0P6 

.600-006 

2.060*010 

25.83 

20 . 0 

-. 100-005 

.600-006 

2.442*010 

30.49 

10.0 

-. 100-006 

. 100-005 

2.121*010 

26.48 

20.0 

-. 100-005 

. 100-005 

2.661*010 

33.22 

IP  . P 

-.2 36-006 

. 100-006 

1 .009*01  0 

23.58 

40  . P 

-. 1 00-006 

. 1 00-006 

3 . 191*009 

3 . 98 

IP  . P 

-.2  36-006 

.236-006 

2. 346*010 

29.28 

40 . 0 

-. 100-006 

. 2 36-006 

4 . 329*009 

5.41 

10.0 

-.236-006 

.600-106 

2.629*010 

32.83 

40  . 0 

-. 100-006 

. 600-006 

5. 148*009 

6.41 

10.0 

-.236-006 

. 100-005 

2.714*010 

33.09 

40.0 

- . 1 00-006 

. 1 00-005 

5.4 14*009 

6 . 76 

10.0 

-.600-006 

. 100-006 

2.060*010 

25.83 

40 . 0 

-.  2 36-006 

. 100-006 

4 . 329*009 

5.41 

10.0 

-.600-#ib 

.236-006 

2.629*010 

32.83 

40 . 0 

-.236-006 

. 2 36-006 

6 . 729*009 

8 .40 

10.0 

-.600-006 

. 600-006 

2.991*010 

37.  14 

40 . 0 

-.236-006 

. 600-006 

8 . 939*009 

11.16 

10.0 

-.6P0-0P6 

. 100-005 

3. 101*010 

30.72 

40.0 

-.236-006 

. 1P0-005 

9.772*009 

12.20 

10.0 

100-0#*) 

. 100-006 

2.121*010 

26.48 

40. 0 

-.600-006 

. 1 00-006 

5 . 1 48*009 

6.43 

10.0 

-. 100-005 

. 2 36-006 

2.714*010 

33.89 

40 . 0 

-.600-006 

.236-006 

0.939*009 

11.16 

10.0 

-. 100-005 

. 600-006 

3. 101*010 

38.72 

40 . 0 

-.600-006 

. 600-006 

1.331*010 

16.62 

10.0 

-. 100-005 

. 100-005 

3.221*010 

40.21 

40 . 0 

-.600-006 

. 1 00-005 

1 . 525*0 1 0 

19.03 

20.0 

-. 100-006 

. 1 00-006 

7.776*009 

9.71 

40 . 0 

-. 100-005 

. 1 00-806 

5. 414*009 

6 . 76 

20 . 0 

-. 1 00-P06 

. 236-006 

1 .005*010 

12.55 

40 . 0 

-. 100-005 

.236-006 

9.772*009 

12.20 

20 . 0 

-. 100-006 

.600-006 

4 • 

40 . 0 

-. 100-005 

.600-006 

1 . 525*01  0 

19.0  1 

20  . P 

-. 100-006 

. 100-005 

1 . 204*010 

15.0  3 

40 . 0 

-. 100-005 

. 1 00-005 

1 . 784*010 

22.27 

20.0 

-.236-006 

. 1 00-006 

1 .005*010 

12.55 

20.0 

-.  2 36-006 

.236-006 

1 .422*010 

17.75 

20.0 

-.236-006 

. 600-006 

1 . 745*010 

21.70 

20  . 0 

-.236-006 

.100-005 

1.053*010 

23.14 

20.0 

-.600-006 

. 100-006 

1 .157*010 

14.44 

20.0 

-.600-006 

.236-006 

1 . 745*010 

21  . 78 

20.0 

-.600-006 

. 600-006 

2 . 2 56*0  1 0 

28.17 

20 . 0 

-.600-006 

. 100-005 

2.442*010 

30.49 

i \bi  i 3ih» 


1 7()U 

al 

236  X I0_,>  m 1 

Ki  ll/ 1 

15  ■ 

Itr1  R.ul 

hull/ 1 

50  / 10"  ' K.ul 

/| 

2oo  in 

, - 

M, 

ni  i m 

»*</>•  «// 

10.  * 

1 00-006 

. 100-006 

6 . 486*009 

27.04 

20 . 0 

-. 1 00-006 

. 100  006 

1 .081*010 

4 4.98 

i0.i) 

-. 1 00-006 

. 236-006 

H. 666*009 

16.61 

20.0 

-. 1 00-006 

. 2 16  006 

1 . 744*010 

72.69 

10.0 

100-006 

.600-006 

9.966*009 

4 1.47 

20 . 0 

-. 100-006 

. 600  -00b 

2.40^*010 

99.97 

10.0 

- . 1 00-006 

. 100-006 

l .041*010 

4 ) . 32 

20 . 0 

-. 100-006 

100-00S 

2.641*010 

109.90 

10.0 

- . 2 16-006 

. 100-006 

8 . 666*009 

36.61 

40 . 0 

-. 100-006 

. 1 00-006 

3. 1 77*009 

l 3.22 

10.0 

-.2)6-006 

.236-006 

i . 263*0  1 0 

62.13 

40 . 0 

-. 100-006 

. 2 36-006 

4 109*008 

17.93 

10  0 

-.2)6-006 

. 600-006 

1.680*010 

66.  74 

40 . 0 

-. 100-006 

6 00  -006 

6.124*00) 

21.3) 

10.0 

-.2)6-006 

. 100-006 

1 .694  *0  10 

70.62 

40 . 0 

-. 100-006 

. 10  0 — 00*' 

6.  188*008 

. . 4. 

10.0 

- .600-006 

. 1 00-006 

9.966*009 

41  . 47 

40.0 

>.  2 36-006 

. 100-006 

4 . 309*008 

17.93 

10.0 

-.600-006 

. 2 36-006 

1 .680*010 

66.74 

40 . 0 

-. 236-006 

. 2 36-006 

6 .697*00') 

27.87 

10 . 0 

-.600-006 

. 600-006 

2. 1 38*010 

88.97 

40 . 0 

-.  2 36-006 

.600-006 

8 .896*009 

37.02 

10.0 

- .600-006 

. 100-006 

2.  364*010 

97.96 

40.0 

-.2)6-006 

. 100-006 

9. 722*009 

40.46 

10.0 

- . 100-006 

. 100-006 

1.041*010 

43.32 

40  . e 

- . 600-00b 

. 100-006 

6 . 1 24*009 

21.33 

10.0 

-. 100-406 

. 2 36-006 

1 .694*010 

7J.  62 

40.0 

- . 600-00b 

. 2 36-006 

8 .896*009 

37.02 

10.0 

- . 1 00-406 

.600-006 

2. 364*010 

97.96 

40 . 0 

- . b00-00b 

. 600-006 

1 . 324*010 

56.11 

10.0 

-.  1 00-00') 

. 100-006 

2.618*010 

108.94 

40 . 0 

-.600-006 

. 100-006 

1 .616*010 

6 3.10 

*0 . 0 

-. 100-006 

. 100-006 

6 . 780*009 

28.22 

40 . 0 

-. 100-006 

. 100-806 

6. 388*009 

22.4  2 

2 0 . 0 

-. 100-006 

. 236-006 

8.906*009 

37.06 

40 . 0 

-.108-006 

.236-006 

9. 722*009 

40 . 4b 

.'0. 0 

-. 1 00-006 

. 600-006 

1 . 0)6*010 

43  . 08 

40.0 

-. 100-006 

. 600-006 

1 .616*010 

6 1.10 

20  . 0 

-. 100-006 

. 100-006 

1 .081*010 

44.98 

40 . 0 

-. 100-006 

. 1 00-006 

1 . 7 74  *0  1 0 

7 3.81 

20  . 0 

-.2)6-006 

. 100-006 

8 . 906*009 

37 . 06 

20.0 

-.2)6-006 

. 236-006 

1 . 297*010 

6 3.98 

2$  . 0 

- . 2 )6-006 

. 600-006 

1 .629*010 

67.78 

2 0 . 0 

-.  2 36-006 

. 100-006 

1.744*010 

72 . 69 

t 0 . 0 

-.600-006 

. 1 00-006 

1 . 0 16*010 

43.08 

20. 0 

-.600-006 

. 2 36-006 

1 .629*010 

67.78 

4 0.0 

-.600-00b 

. 600-006 

2. 188*010 

91  . 06 

00.0 

- . 600-006 

. 100-006 

2. 402*010 

99 .97 

32 
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l = 700.  a t = 23(1  X I0“6  m_l 


Re 

(Wl  = 1 5 X 

10*1  R.td 

Unit) ) 

- 15  X 

l(T:  K iel 

Z,  = 200. 

111 

0lv»* 

(HAl  M„* 

88 'i- 

Of  l M 

1 HAi  m / 

re  • 

f 

-.  1 00-006 

. 100-006 

2 . 172*009 

27.12 

29. 9 

- . 100-005 

. 100-006 

3.614*009 

45.12 

li.  8 

100-006 

.236-006 

2.861  *009 

35.72 

20 . 0 

-.  100-005 

. 2 16-006 

5.833*009 

72.83 

10.0 

100-006 

. 600-006 

3 . 3 32*009 

41.60 

20 . 0 

-.  100-005 

. 600-006 

8 . 0 33*009 

180.29 

10.0 

100-006 

. 10  0 -00' 

3.401*009 

4 3.46 

20 . 0 

-.  1 00-005 

. 100-005 

8.904*009 

111.17 

10.0 

- . 2 36-006 

. 100-006 

2.861*009 

35.72 

40.0 

- . 100-006 

. 100-006 

2 . 506*009 

11.28 

10.0 

-.236-006 

.236-006 

4 . 1 8 9 *0 0 9 

52.  30 

40 . 0 

-.  100-006 

.236-006 

3.268*009 

40.81 

i0.0 

-.236-006 

. 600-006 

5 . 28 2 *009 

65.95 

40.0 

-.  100-006 

. 600-006 

3. 702*009 

47.22 

10.0 

-.236-006 

. 100-005 

5.666*009 

70  . 74 

40 . 0 

- . 1 0 0-006 

. 100-805 

1.943*009 

49.23 

10.0 

-.600-006 

. 100-006 

3 . 3 32*009 

41  .60 

40 . 0 

-.2  36-006 

. 100-006 

3 . 268*009 

40.81 

10.0 

-.600-006 

.236-006 

5.202*009 

65.95 

40.0 

-.236-006 

.236-006 

4 . 700*009 

58.67 

10.0 

- . 600-006 

.600-006 

7.  1 49*009 

89 . 26 

40 . 0 

-.236-006 

.600-806 

5.840*009 

72.91 

10.0 

-. 600-006 

. 100-085 

7 .870*009 

98 . 26 

4 0.0 

- . 2 36-00C 

. 1 00-085 

6.212*009 

77.81 

10.0 

-.  100-805 

. 1 00  - 0 0 b 

3 .481  *00  9 

4 3.46 

40.0 

-.600-006 

. 100-006 

3. 782*009 

47.22 

10.0 

100-005 

. 2 36-006 

5.666*009 

70.  74 

40.  0 

-.600*006 

.236-006 

5 . 840*009 

72.91 

10.0 

- . 100-005 

. 600-006 

7.870*009 

98 . 26 

40.0 

-.600-006 

. 600-006 

7. 712*009 

96 . 20 

10.0 

-. 100-005 

. 100-005 

8.  754*009 

109.29 

40 . 0 

- . 600-006 

. 100-005 

8.411*009 

185.01 

29. i 

-. 100-006 

. 100-006 

2.267*009 

28.31 

40 . 0 

-.  100-805 

. 100-006 

3.94  3*089 

49.23 

29 . 0 

-. 100-006 

. 236-006 

2 . 978*009 

37.18 

40 . 0 

-.100-005 

.236-006 

6 .212*009 

77.81 

20.0 

-.  1 00-006 

. 600-006 

3.462*009 

4 3.22 

40 . 0 

-.  1 00-005 

.600-006 

8.41 1*009 

105.01 

20.0 

-.  1 00-006 

. 100-00*) 

3.614*009 

45.  12 

40 . 0 

-.100-005 

. 100-005 

9.249*009 

115.48 

20.  0 

-. 2 36-006 

. 100-006 

2 .978*009 

37.18 

29.9 

-.236-006 

.236-006 

4 . 337*009 

54.15 

29.9 

-.2  36-006 

. 600-006 

5 .446*009 

68 . 00 

29.4 

-.2  36-006 

. 100-005 

5.833*009 

72.83 

29  9 

-.600-006 

. 100-006 

3 . 462*009 

41.22 

20  0 

-.600-006 

. 2 36-006 

5.446*009 

68.00 

20  0 

-.600-006 

.600-006 

7. 317*009 

91  . 36 

i0  0 

-.600-006 

100-00*) 

6 .0  3 1*009 

100.29 

i \m  i Kii) 


1 '(III 

a j = 230  X 1 0 ’ m ~ 

Ri  in  i 

15  ■ 

' l(T  1 K.id 

ImOM  50  / 1 0“ - R.kJ 

:on  m 

i.. . 

n»  i ** 

«,./  H* 

10.0 

- . i 99-006 

. : #0-006 

6 . 767*008 

28.12 

20.0 

- . 2 00-006 

. 1 00-006 

1.126*009 

4 o . 86 

* 

- . 100-006 

. 2 16-006 

8.901 *038 

17  . 04 

20.0 

- . 1 00-006 

. 2 36-006 

1 .817*009 

76.62 

it.# 

. 1 #0-006 

. 600-006 

1 . 0 3 7 *099 

4 3.16 

20.0 

-. 100-006 

. 600-006 

2 . 601*009 

104.16 

:# . # 

- . 1 00  -006 

. 1 00-00S 

1 .081*009 

46.07 

20 . 0 

-. 100-006 

. 1 00-006 

2 . 774*009 

116  46 

\i.  0 

-.2  16-006 

. 1 00-096 

8. 9014008 

04 

40 . 0 

-.  00-006 

. 100-006 

7.818*008 

12 . 64 

it.  # 

-.2  16-006 

. 2 36-006 

1 . 301*099 

4 . 

40.0 

-. 1 00-006 

. 2 36-006 

1 .020*009 

42.44 

l# . 0 

-.216-906 

.600-906 

1 .044*009 

68.42 

40 . 0 

-. 100-006 

.600-006 

1 . 1 80*009 

49.11 

li.# 

- . 2 36-006 

. 1 00 -00  S 

1 . 764*909 

7 1.39 

40 . 0 

-. 1 00-006 

. 100-086 

1 . 2 30*009 

61  .19 

i#.# 

- . 6 #0-006 

. 1 00-096 

1 .017*009 

4 3.16 

4?  . 0 

-.2  16-006 

. 1 00-006 

1 .020*009 

42  . 44 

it  # 

- . 609-006 

. 2 36  — 006 

1 . 644*029 

68.42 

40.0 

-.236-006 

.236-006 

1 . 466*009 

61.01 

l#.  # 

- .000-006 

. 690-006 

2.226*009 

92.6) 

40 . 0 

-.2 16-006 

. 600-006 

1 .821  *009 

76.81 

it.# 

-.600-006 

. 100-006 

2.460*009 

102 . 98 

40.0 

-.  2 36-006 

. 200-006 

1 . 944  *009 

80.89 

li.e 

-.  109-00S 

. 100-006 

1 . 08  1*009 

46.07 

40 . 0 

- . 600-006 

. 1 00-006 

1 . 180*009 

49.11 

19.  # 

- . 100-00S 

. 2 36-006 

1 . 764  *009 

7 1.39 

40 . 0 

-.600-006 

.216-006 

1.821 *009 

76.81 

1#.# 

-.  l#0-0#6 

. 000-006 

2.460*009 

101.98 

40 . 0 

- . 600-006 

. 600-006 

2.406*009 

100.08 

it.# 

- . 1 00-006 

. 1 00-00* 

2. 726*009 

111.44 

40 . 0 

-.600-006 

. 1 00-006 

2.622*009 

109.14 

it . 9 

- . i 00-006 

. 100-006 

7.062*008 

29.  19 

40 . 0 

-. 100-006 

. 1 00-006 

1 . 2 10*009 

61.19 

29.  e 

100-006 

. 2 16-006 

9 . 276*008 

18 . 60 

40 . 0 

-. 100-006 

. 2 36-006 

1 .944*009 

80 . 89 

it  . # 

1 00-006 

. 6 0 0-006 

1 .078*009 

44.88 

40 . 0 

-. 100-006 

. b00 -006 

2.622*009 

109.14 

it  ■ 9 

100-006 

. 100-906 

1 . 1 26*009 

46 . 86 

40 . 0 

-. 100-006 

. 1 00-006 

2 .884*009 

120.01 

2t . 9 

-.216-096 

. 100-006 

9 . 276*008 

18.60 

.9  9 

-.2  16-006 

. 2 36-006 

1 . 1614009 

66.21 

2# . 0 

-.216-006 

. 600-096 

1 .69^*009 

70.61 

<’0.0 

-.236-006 

. 100-006 

1.81  7*009 

76.62 

<0.0 

-.600-006 

. 1 00-006 

1 .070*009 

44.88 

29.0 

-.600-006 

. 2 16-006 

l .697*009 

70.61 

it  . 9 

-.600-006 

. 600-006 

2.280*009 

94.87 

;t . 0 

-.600-006 

. 100-006 

2.60  1*009 

104.16 

i4 


TABU  5<o> 


1 = 700  a,  = 236  X I0_(>  m'1 


Kc 

( ) = .|S  X 

10  ■ Rad 

lm(0 ) 

= 15 

X 10  • Rad 

/ ! = 500 

111 

i i <*>  ) 

...  . 

Dfi  » 

t M/,  , M«, 

AIM, 

10.0 

- . 1 30-306 

. 1 00-006 

6 .929*009 

8 . 66 

20.0 

-. 100-306 

. 100-306 

6 .092*009 

7.61 

10.0 

-.  1 00 -00b 

. 2 36-006 

7 . 4 6 2 ♦ 0 3 9 

9.32 

20.0 

-.100-006 

.236-006 

7 . 406*039 

9.25 

10.0 

-.  100-006 

. 600-006 

7.  727*009 

9.65 

20.0 

-.100-006 

. 600-006 

8. 1 95*039 

10.23 

10.0 

-. 100-006 

. 100-006 

7 .799*009 

9.74 

20.0 

-.103-305 

. 100-006 

8.428*009 

10.62 

10.0 

- . 2 36-006 

. 100-306 

7 . 462*009 

9.  32 

40 . 0 

-. 100-006 

. 100-306 

2.535*009 

3.17 

10.0 

-. 236-006 

.236-006 

6.064*009 

10.09 

40 . 0 

-. 103-006 

.236-006 

3.204*009 

4.00 

10.0 

-.236-006 

. 600-006 

8.  396*009 

10.48 

40.0 

-.100-006 

. 600-006 

3.632*009 

4 .53 

10.0 

-.236-006 

. 100-006 

8 . 481  *009 

10.59 

40 . 0 

-.100-006 

.103-006 

3 . 762*009 

4 . 70 

10.0 

-.600-006 

. 100-006 

7. 727*009 

9.66 

40.0 

-.236-006 

. 1 00-006 

3.204*009 

4 . 00 

10.0 

-.600-006 

.236-006 

S. 396*009 

10.48 

40 . 0 

-.  2 36-006 

.236-006 

4 . 354*009 

6.44 

10.0 

-.600-006 

. 600-006 

8 . 73  3*009 

10.90 

40.0 

-.  2 36-036 

.600-306 

5.183*309 

6.47 

10.0 

-.600-006 

. 100-006 

8 . 826*009 

11.02 

40.0 

- .236-006 

. J 03-005 

5.452*009 

6.81 

10.0 

-. 100-006 

. 100-006 

7.799*009 

9.74 

40.0 

- .600-336 

. 100-306 

3.632*009 

4.53 

10.0 

-.  100-006 

.236-006 

8.481*009 

10.59 

40 . 0 

-.600-006 

.236-036 

6. 183*009 

6 .47 

10.0 

- . 100-006 

.600-006 

8.826*009 

11.02 

40 . 0 

- .600-006 

. 600-006 

6. 402*309 

7.99 

10.0 

-.  100-006 

. 100-006 

8.919*009 

11.14 

40.0 

-.630-006 

. 100-336 

6 .818*089 

8.51 

20.0 

-.  1 00-006 

. 100-006 

4 . 769*009 

5.95 

40 . 0 

-.  1 0 0-00-,, 

. 1 00-816 

3 . 762*039 

4 . 70 

20 . 0 

-.100-006 

.236-006 

6.639*009 

6.92 

40 . 0 

-. 100-036 

.236-306 

5.452*039 

6 .61 

20.0 

-.  100-006 

.600-006 

6.969*009 

7.45 

40 . 0 

-.133-336 

. 600-306 

6 . 818*009 

6.51 

20 . 0 

-.  100-006 

.100-006 

6.092*009 

7.61 

40.0 

- .103-006 

. 100-305 

7.292*009 

9.10 

20 . 0 

-.2  36-006 

. 100-006 

6.639*009 

6.92 

20 . 0 

-.236-006 

.236-006 

6 . 606*009 

8.25 

20.0 

-.236-006 

. 600-006 

7.226*089 

9.02 

20 . 0 

-.236-006 

. 100-006 

7 . 406*009 

9.25 

20.0 

-.600-006 

. 100-006 

6.969*009 

7.45 

20.0 

-.600-006 

.236-006 

7.226*009 

9.02 

20 . 0 

-.600-006 

.600-006 

7 .976*009 

9.96 

20.  0 

-.600-006 

. 100-016 

8.  196*009 

10.23 

I AH  I I id i 


1 700 

a j 2. ill  / 1 (1  (l  in  ^ 

Hcu)  i 

15  ■ 

■ 10*  1 K .nl 

ImlH  I 50  X 10“  ' K.iJ 

/.  | 500  in 

• “ 

t>*  i M 

10.0 

- . 1 00-006 

. 1 00-006 

6 . 746*009 

28.0  7 

20.0 

-. 100-005 

. 1 00-006 

6.062*009 

25.23 

10  .0 

- . 1 f0  - 0 0 fc 

. 2 36-006 

7 . 420*009 

30 . 88 

20 . 0 

-. 100-005 

. 2 36-006 

7.377*009 

30  . 70 

10.0 

-.  1 00-006 

.600-006 

7.687*009 

31.99 

20.0 

-. 1 00-005 

.600-006 

8. 166*009 

3 3.99 

.0.1* 

-.  100-006 

. 100-005 

7 . 759*009 

32.29 

20.0 

-. 100-005 

. 100-005 

6.  399*009 

34 . 96 

10.0 

-.2  36-006 

. 100-006 

7.420*009 

30 . 88 

40 . 0 

-. 100-006 

. 1 00-006 

2 . 5 2 5 * 13  0 9 

10.51 

10  0 

-.236-006 

. 2 36-006 

8 .046*009 

33.49 

40.0 

-. 1 00-006 

. 2 36-006 

3.192*009 

13.29 

10.0 

-.236-006 

. 600-006 

8.361*009 

34 . 80 

40.0 

-.  1 00-006 

. 600-006 

3.619*009 

15.06 

10.0 

-.236-006 

. 100-005 

8 .447*009 

35.15 

40.0 

- . 100-0  06 

. 100-005 

3 . 749*009 

15.60 

10.0 

-.600-006 

. 1 00-006 

7.687*009 

31.99 

4 0 . 0 

2 36-006 

. 100-006 

3. 192*009 

1 3.29 

10.0 

- .000-006 

. 2 36-006 

8 . 361 *009 

34.80 

4 0.0 

- . 2 36-006 

. 2 3b  — 006 

4 . 3 38  *009 

18.05 

10.0 

- . 600-006 

.600-006 

8.701 *009 

36.21 

40.0 

-.  2 36-006 

.600-006 

S. 165*009 

21  .49 

10.0 

-.600-006 

. 100-005 

8.794*009 

36.60 

40 . 0 

-.236-006 

. 100-005 

5.4  3 3 *009 

22.61 

10.0 

-. 1 00-005 

. 100-006 

7 .759*009 

32.  29 

40 . 0 

- . 600-006 

. 1 00-006 

3.619*009 

15.06 

10.0 

1 00-005 

. 2 36-006 

8 . 447*009 

3 c .15 

40.0 

—.600—006 

. 2 36-006 

5. 165*009 

21  .49 

10.0 

-. 1 00-005 

. 600-006 

8.794*009 

36.60 

40.0 

— . 600-006 

. 600-006 

6.381 *009 

26.56 

If  .0 

-. 100-005 

. 100-005 

9 .889*009 

36 . 99 

40 . 0 

- .600-006 

. 100-005 

6 . 796*009 

28.28 

20. 0 

-. 1 00-006 

. 1 00-006 

4 . 74  3 *009 

19.74 

40 . 0 

-.100-005 

. 1 00-006 

3 . 749*009 

15.60 

.0.0 

-. 1 00-006 

. 2 36-006 

5.511*009 

22.93 

4 0.0 

-.100-005 

.236-006 

5.4  3 3 *009 

22.61 

20. 0 

-.  1 0 0 — 0 0 b 

.600-006 

5.940*009 

24  . 72 

40 . 0 

100-005 

. 600-006 

6 . 796*009 

28. 28 

20.0 

-.100-006 

. 1 00-005 

6.062*009 

25.2  3 

40 . 0 

-. 100-005 

. 100-005 

7 . 269*039 

30.25 

20.0 

-.  2 36-006 

. 1 00-006 

5.51 1*009 

22.93 

20 . 0 

-.  2 36-006 

. 2 36-006 

6 .576*009 

27.  37 

20.0 

-.2 36-006 

.600-006 

7. 196*009 

29.95 

20.0 

- . 2 36-006 

. 1 00-005 

7. 377*009 

30. 70 

20. 0 

-.600-006 

. 100-006 

5.940*009 

24.72 

*0.0 

-.600-006 

. 2 36-006 

7. 196*009 

29.95 

20.0 

-.600-006 

. 600-006 

7.946*009 

33.07 

20.0 

- .600-006 

. 100-005 

H . 166*009 

33.99 

I AH!  I ?<gi 


1 = 

700. 

a 

1 - 2M>  X I O''1  m;~ 

1 

Rcltf  ! - 1 5 X 

ur1  k.ki 

II 

\r\ 

X 

Hr2  Rad 

/ , 500. 

rn 

f«T  y. 

ATNU 

:(  ** 

IHAl  M/ 

Lift/*-  ) 

• ** 

Hi 

18/*- t 

10.0 

100-006 

. 100-006 

2.256*009 

28.  16 

20.0 

-.  100-005 

. 1 00-006 

3.861*009 

48.08 

10.0 

1 00-006 

. 2 36-006 

3.007*009 

37.54 

20.0 

-.  1 00-005 

. 2 36-006 

6.475*009 

80 . 84 

10.0 

-.  100-006 

. 600-006 

3.532*009 

44.10 

20 . a 

-.  100-005 

. 600-006 

7.921*009 

98 . 89 

10.0 

1 00  - 0 0b 

. 100-005 

3. 700*009 

46.19 

20 . 0 

-.100-005 

. 1 00-005 

B . 1 57*009 

101.84 

10.0 

236-006 

. 100-006 

3.007*009 

37.54 

40 . 0 

-.  100-006 

. 1 00-006 

2.442*009 

10.49 

10.0 

-.2 )b -0  0 b 

.236-006 

4 .510*009 

5b.  30 

40 . 0 

-.  100-006 

. 2 36-006 

3.089*009 

18.57 

10.0 

-.2 36-006 

. 600-006 

5.601*009 

72.46 

40 . 0 

- . 1 00-006 

.600-006 

3.504*009 

4 3.  74 

10.0 

-. 236-00b 

. 1 00-005 

6 .270*009 

78.28 

40 . 0 

1 00-006 

. 100-006 

3.6)0*009 

45.32 

10.0 

- . 600 -006 

. 100-006 

3.5)2*009 

44 .10 

40 . 0 

-.  2 36-006 

. 100-006 

3 .089*009 

10.67 

10.0 

-.600-00b 

.236-006 

5.601*009 

72.46 

40.0 

-.  236-006 

.236-006 

4 . 204*009 

52.  49 

10.0 

- . 600 -006 

. 600-006 

6 . 1 10*009 

101.60 

40 . 0 

-.  2 36-006 

.600-006 

5.01 1*009 

62.56 

10.0 

-.600-006 

. 100-005 

0.531*009 

106.51 

40 . 0 

-.236-006 

. 1 00-005 

5.273*009 

66.03 

10.0 

100-005 

. 1 00-006 

3 . 700*009 

46.19 

40 . 0 

-.600-006 

. 1 00-006 

3.604*009 

4 1.74 

10.0 

- . 1 00-005 

. 236-006 

6.270*009 

78.28 

40.0 

- .600-006 

. 2 16-006 

5.011*009 

62.56 

10.0 

-. 100-005 

. 600-006 

8.531*009 

106.51 

40 . 0 

- . 600-006 

.600-006 

6 . 200*009 

77.41 

10.0 

-.100-005 

. 100-005 

8.634*009 

107.80 

40 . 0 

-.600-006 

. 100-006 

6.607*009 

82  . 49 

2 0 . 0 

-. 100-006 

. 100-006 

2.358*009 

29.44 

40 . 0 

-.  100-006 

. 1 00-006 

3 . 6 30*009 

4 5.  32 

20 . 0 

100-006 

.236-006 

3. 137*009 

39.  16 

40 . 0 

-.  100-006 

. 2 36-006 

5.273*009 

65.83 

20.0 

-.100-006 

. 600-006 

3.678*009 

45.92 

40.0 

-.  1 00-005 

.600-006 

6.607*009 

82  . 49 

20.0 

-.100-006 

. 100-005 

3.851*009 

48.08 

40 . 0 

- . 1 00-006 

. 1 00-005 

7.071*089 

08 . 28 

20.0 

-.  2 36-006 

. 100-006 

3. 137*009 

39. 16 

20 . 0 

-.236-006 

.2)6-006 

4 .602*009 

58 . 46 

20.0 

-.236-006 

. 600-006 

6 . 002*009 

74.93 

20 . 0 

-.236-016 

. 100-005 

6 .475*009 

80 . 04 

20.0 

-.600-006 

. 1 00-006 

3.678*009 

45.92 

20 . 0 

-.600-006 

. 2 36-006 

6 . 002*009 

74.93 

20 . 0 

-.600-006 

. 600-006 

7.698*009 

96.11 

20 . 0 

-.600-006 

. 100-005 

7.921*009 

98 . 89 

37 


I \M  I ?<ln 


1 7(lt) 

f,l 

X Kr'1  m 1 

Ki'i  (h 

15  ' 

nr1  K.ni 

lmt(i  i 

su  x nr-  k.hI 

/i 

500.  m 

“ 

m- 

« 

‘ » * .• 

i e . 9 

- . [00-00* 

. 100-006 

7.016*008 

29. 20 

29.0 

-. 100-206 

. 1 00-006 

1 . 1 99*009 

49.91 

10.0 

- . 2 00-00b 

. 2 16-006 

9.  164*00* 

18.93 

29.9 

-. 1 00-006 

. 2 36-006 

2.01 7*009 

8 3.94 

10  . 0 

-. 100-006 

.600-006 

1 .0*49*00* 

46.7) 

2 9.9 

- . 100-006 

.600-006 

2.898*009 

120.60 

19.9 

-. ; 00-006 

. 1 00-00*' 

1 . 161*009 

4 7 . 90 

29 . 9 

- . 100-006 

. 1 00-006 

3 . 268*009 

1 36 . 01 

10.0 

- . 2 l*-00* 

. 100-006 

9.  364*000 

3 8.91 

40.0 

- . 1 0 0—00** 

. 100-006 

8. 166*008 

33.98 

19.9 

- . 2 36-006 

. 2 36-006 

1 . 40  3 *009 

68.  19 

40 . 0 

100-006 

. 2 36-006 

1 .080*009 

44.94 

10.0 

- . 2 36-006 

. 6 0 0-006 

1 . 806*009 

76.16 

40  0 

-. 1 00-006 

.600-006 

1 . 261 *009 

62.48 

10.9 

-.2  lb-09* 

. 100-006 

1 . 961*009 

81.20 

40.0 

- . 1 00-006 

. 100-006 

1 . 318*009 

64.87 

1 9.9 

- . *99-09* 

. 100-006 

l .099*009 

4 6.73 

40 . 0 

-.  2 36-006 

. 100-006 

1 .080*009 

44.94 

10.0 

-.*00-00* 

.2 36-00b 

1 . 806*009 

76.16 

40 . 0 

-.236-006 

. 2 16-006 

1 . 693*009 

66.31 

10.0 

- .600-006 

. 600-006 

2.6  3 3 *009 

106.41 

40 . 0 

-.  2 36-006 

.600-006 

2.022*009 

84  . 16 

19.0 

-.000-006 

. 100-007 

2.828*009 

117.70 

40 . 0 

-.  2 36-006 

. 100-006 

2. 1 74*009 

90 . 48 

10.0 

- . 1 00-00*) 

. 100-006 

1 . 161*009 

47.90 

40 . 0 

- .600-006 

. 1 00-006 

1 . 261*009 

62 .48 

20.0 

-.  1 00-006 

. 2 36-006 

l .961  *009 

81.20 

40 . 0 

-.600-006 

. 2 36-006 

2.072*009 

84.16 

if  .0 

-.  [00-90,> 

. 600- 006 

2.828*00 9 

117.70 

40. 0 

-.600-006 

.600-006 

2. 767*009 

116.17 

10 .0 

-.  100-006 

. 100-006 

3.201 *009 

113.2  3 

40 . 0 

-.600-006 

. 1 00-006 

3.060*009 

127.34 

*9  0 

-. 100-006 

. 1 00-006 

7 . 346*008 

30.67 

40 . 0 

-. 100-006 

. 1 00-006 

1 . 318*009 

64.07 

* 0 . * 

- . 100-006 

. 2 36-006 

9. 769*008 

40 . 66 

40 . 0 

-. 100-006 

. 2 36-006 

2.174*309 

90 . 48 

*0 . 0 

-. 100-006 

. 600-006 

1 . 146*009 

47.68 

40.0 

-.  100-006 

. 600-006 

3 . 060  *009 

127.34 

10 . 0 

- . i 00-006 

. 1 00-006 

1 . 199*009 

49.91 

40.0 

-. 100-006 

. 100-006 

3.421*009 

142.  38 

<0 . 0 

- . 2 16-006 

. 1 00-006 

9 . 769*008 

40.66 

20 . 0 

- . 2 16-006 

. 2 36-006 

1 .468*009 

60.70 

20.0 

-. 236-006 

. 600-006 

1 . 869*009 

77.80 

20 . 0 

- . 2 16-006 

. 100-006 

2.01 7*00* 

83.94 

29 . 9 

-.600-006 

. 100-006 

1 . 1 46*009 

47.68 

29 . 9 

- .600-006 

. 2 36-006 

1 . 869*009 

77 . 80 

20 . 0 

-.600-006 

. 600-006 

2.603*009 

108.33 

20 . 0 

- . 600-006 

. 100-006 

2 . 898*009 

120.60 

.IS 


TABU  l(n 


I = 700  Q | = :.0>  X 10  111-1 


Ri'UM  = 1 5 X 

10 K.iil 

Imli/ ) 

- 15  ■ 

10~*  R.id 

i 

f | 1000 

111 

fHA*  •*.- 

' • 

Of  . *• 

, MAX  H/ 

AT*  i 

ie,. 

19.0 

-. 1 00-006 

. 1 00-006 

2. 910*009 

3.64 

20.0 

100-006 

. 1 00-006 

2.760*009 

3.44 

If  .0 

100-006 

.236-006 

3 . 000*009 

3.  76 

20.0 

-. 100-005 

.236-006 

2.999*009 

3.  74 

If  . 0 

-.100-006 

.600-006 

3.050*009 

3.81 

20.0 

-. 100-005 

.600-006 

3. 121*089 

3. 90 

lf.0 

-. 100-006 

. 100-005 

3.061*009 

3.82 

20 . 0 

-. 100-005 

. 100-006 

3. 164*009 

3.94 

10.0 

-.236-006 

. 100-006 

3 .008*009 

3.76 

40.0 

-. 100-006 

. 100-006 

1 .687*009 

2.1  1 

10.0 

-.  2 36-006 

. 2 36-006 

3 . 1 04  *009 

3 . 88 

40 . 0 

100-006 

2 36-006 

1 .969*009 

2.46 

10.0 

-.2  36-006 

. 600-006 

3. 149*009 

3.91 

4 0.0 

-. 100-006 

. 600-006 

2.1 11*809 

2.64 

If  . 0 

-.236-006 

. 100-006 

3.161 *009 

3. 95 

40 . 0 

100-006 

. 100-005 

2 . 1 64*009 

2.69 

10.0 

-.600-006 

. 100-006 

3 .050*009 

3.81 

40 . 0 

-.  2 36-006 

. 100-006 

1 . 969*009 

2.45 

10.0 

-.600-006 

.236-006 

3 . 1 49*009 

3.93 

40.0 

-.  2 36-006 

. 236-006 

2. 336*009 

2.92 

10.0 

-.600-006 

.600-006 

3. 196*009 

3 . 99 

40.0 

-.  2 36-006 

.600-006 

2.566*009 

3.19 

10.0 

-.600-006 

. 100-005 

3 .208*009 

4 . 00 

40.0 

-.236-006 

. 100-006 

2.619*009 

3.27 

10.0 

-. 100-005 

. 1 00-006 

3.061*009 

3.02 

40 . 0 

-.600-906 

. 100-806 

2.11 1*009 

2.64 

10.0 

-. 100-005 

.236-006 

3.  161*009 

3.96 

40.0 

-.600-006 

. 2 36-806 

2.665*009 

3.19 

20.0 

- . i 00-005 

.600-006 

3. 208*009 

4 . 00 

40.0 

-.600-006 

. 600-006 

2.020*009 

3.52 

10.0 

-. 1 00-005 

.100-005 

3.220*009 

4.02 

40 . 0 

-.600-006 

. 100-006 

2 . 090*809 

3.62 

2f.  0 

-. 100-006 

. 100-006 

2.450*009 

3.06 

40 . 0 

-.100-005 

. 100-006 

2. 1 54*089 

2.69 

20.0 

-. 100-006 

.236-006 

2 . 6 39*009 

3.29 

40 . 0 

-.100-005 

.236-006 

2.619*009 

3.27 

20.0 

-. 100-006 

.600-006 

2.733*009 

3.41 

40 . 0 

100-005 

.600-006 

2 .890*009 

3.62 

20. 0 

-. 100-006 

. 100-005 

2.758*009 

3.44 

40 . 0 

-. 100-006 

. 100-006 

2.900*809 

3.72 

20.0 

-.236-006 

. 100-006 

2.639*009 

3.29 

20.9 

-.236-006 

.236-006 

2. 859*009 

3.57 

20 . 0 

-.236-006 

. 600-006 

2.969*009 

3.71 

20 . 0 

-.236-006 

. 100-005 

2 . 999*009 

3.74 

20.0 

-.600-006 

. 1 00-006 

2.733*009 

3.41 

20.0 

-.600-006 

.236-006 

2 . 969*009 

3.71 

20.0 

-.600-006 

. 600-006 

3.088*009 

3.86 

20.0 

-.600-006 

. 100-005 

3.121*009 

3.90 

I Mil  I ii|i 


1 = 700 

° i 

: 5(.  x |()-(>  m 1 

RilR  l 

15  > 

nr1  k.k! 

fm  < i 

50  X 10'  ' R.hI 

/l 

loot)  111 

PI  - 

»*••*  • 

,l  i ** 

"• 

if* . 0 

-. 1 00-006 

. 100-006 

2 . 909*004 

12.10 

20.0 

-.  100-005 

. 10P-006 

2. 751*009 

1 l . 45 

l* .0 

100-006 

. 2 36-006 

3 . 080*009 

12.49 

20.0 

-. 100-005 

. 2 36-306 

2.991*009 

12.46 

10.  0 

100-00b 

. 600-006 

3 .04 3*009 

12.66 

20.0 

-. 100-005 

. 600-006 

3. 1 15*009 

12.96 

10.0 

- . 1 P0-00b 

. 100-085 

3 . 0 5 4 ♦ 0 0 9 

12.71 

20.0 

-. 100-805 

. 1 00-005 

3. 149*009 

13.18 

10.0 

-.2  36-006 

. 1 00-086 

3 . 000*009 

12.49 

40.0 

-. 100-006 

. 100-006 

1.682*009 

7 . 00 

10.0 

-.2  36-086 

. 2 36-006 

3 . 098*009 

1 2.H9 

4 0 . 0 

-. 100-006 

. 2 36-006 

1 .954*009 

8.13 

10.0 

-.2  36-006 

. 600-006 

3. i 43*009 

13.08 

40.0 

100-006 

. 600-006 

2. 106*009 

0 . 76 

10.0 

-.23 6-006 

. 100-00*-. 

J. 1 55*009 

13.13 

40 . 0 

-. 1 00-006 

. 1 00-305 

2. 1 49*009 

8 .94 

10.0 

- . 600-006 

. 1 00-006 

3 .04  3 *009 

12.66 

40 . 0 

-.  2 36-006 

. 100-006 

1.954*009 

0.13 

10.0 

-.b00-00b 

. 236-006 

3. 1 43*009 

13.08 

40 . 0 

-. 236-006 

.236-306 

2. J 3 1 *009 

9.  70 

10.0 

-.600-006 

.600-006 

3. 190*009 

1 3.  28 

40.0 

-.  2 36-006 

. 600-006 

2.550*009 

10.61 

10.0 

- . b00-00b 

. 1 00-005 

J . 202*009 

1 3.  33 

48 . 0 

-.  2 36-006 

. 100-005 

2.614*009 

18.80 

10.0 

-.100-005 

. 1 00-006 

3.854*009 

12.71 

40 . 0 

-.600-006 

. 100-006 

2. 106*009 

8.  76 

10.0 

-. 100-005 

. 2 36-006 

3. 155*009 

13.13 

40 . 0 

-.600-006 

. 2 36-006 

2.550*009 

10.61 

10.0 

-.  100-005 

. 600-806 

3 . 282*009 

13.33 

40.0 

-.600-006 

.600-006 

2 .815*009 

11.71 

10.0 

— . 100-005 

. 100-005 

3.215*009 

1 3.  38 

40 . 0 

-.600-006 

. 100-005 

2.893*089 

12.04 

20. 0 

-.  1 0 0 - 0 0 b 

. 100-006 

2.443*009 

10.17 

40.0 

-. 100-005 

. 100-006 

2. 149*009 

0. 94 

20 . 0 

-. 100-006 

. 2 36-006 

2.6  32*009 

10.95 

40.0 

-. 100-305 

. 2 36-086 

2.614*009 

10.88 

20 . 0 

-. 100-008 

. 600-006 

2 . 726*809 

11.34 

40 . 0 

-. 100-005 

. 600-006 

2.893*009 

12.04 

20 . 0 

1 00 -00b 

. 100-005 

2 . 751*809 

11.45 

40 . 0 

100-005 

. 100-006 

2.975*009 

12.38 

; 0 . 0 

-. 236-006 

. 100-806 

2 . 632*009 

10.95 

2 0 . 0 

-.236-006 

. 2 36-006 

2 . 852*009 

11.87 

20.0 

-.236-006 

. 600-006 

2 . 96  3*009 

12.33 

20.0 

-.  2 36-006 

. 100-005 

2. 99  3*009 

12.46 

20.0 

-.600-006 

. 100-006 

2 . 726*009 

11.34 

20 . 0 

- .600-006 

. 2 36-006 

2.963*009 

12.33 

.-0.0 

- .600-006 

. 600-006 

3.083*009 

12.83 

20.0 

-.600-006 

. 100-005 

3. 1 15*009 

12.96 

40 


I 


I ABU  mi 

I = 700.  a I 23f>  X \0~(' m-1 


K;' 

Hi)  = 1 5 X 

10"1  K.kI 

Imltf  i 

X 

yr. 

II 

l(T:  K.kI 

(.  | 1000 

111 

Mi" 

I 9 0 

1 00-006 

. 100-006 

2.639*009 

ATh 

31.70 

nf  i * 

20 . 0 

-.  1 00-006 

. 1 00-006 

2.691 *009 

3 3.69 

10.0 

-.  100-006 

. 2 36-006 

2.931 *009 

36.60 

20.0 

- . 100-006 

. 2 36-006 

2.941 *009 

36  . 72 

10.0 

- . 1 00-006 

. 600-006 

2. 97# *009 

37.18 

20. 0 

-.  1 00-006 

. 600-006 

3. 060*009 

38.  30 

10.0 

1 00-00b 

. 100-006 

2 . 991*009 

37.  34 

20 . 0 

- . 100-006 

. 100-006 

3.10  1*009 

38  . 74 

10.0 

-.236-006 

. 100-006 

2 .931*909 

36 . 60 

40 . 0 

- . 100-006 

. 1 00-806 

1.641 *009 

20.49 

10.0 

-.2  36-006 

. 2 36-006 

3 . 039*009 

37.94 

40.0 

- . 1 00-006 

. 2 36-006 

1.911 *009 

2 3.06 

10.0 

-.236-006 

. 600-006 

3 . 089*009 

38.67 

40.0 

- . 100-006 

. 600-006 

2.061 *009 

26.  74 

10.0 

-.236-006 

. 100-006 

3. 103*009 

38.  74 

40 . 0 

- . 1 0 0 - 0 0 b 

. ] 00-086 

2. 104*009 

26 . 27 

10.0 

-.600-006 

. 100-006 

2.976*009 

37.18 

40  . 0 

- . 2 36 -00h 

. 100-006 

1.911 *809 

2 3. 06 

10.0 

- .600-006 

.236-006 

3 .089*009 

38.67 

40 . 0 

-.  2 36-806 

. 2 36-006 

2. 286*009 

28.63 

10.0 

- . 600-006 

. 600-006 

3.142*009 

39.22 

40 . 0 

- . 2 36-006 

. 600-806 

2.606*089 

31.27 

10.0 

-.600-006 

. 100-006 

3 . 1 66*009 

39.40 

40 . 0 

-.2  36-006 

. 100-006 

2 . 668*809 

32.07 

10.0 

-. 100-006 

. 1 00 -00b 

2.991 *009 

37.  34 

40 . 0 

- . 600-006 

. 1 00-806 

2.061*009 

26.74 

10.0 

-.  1 00-006 

.236-006 

3 . 1 0 3 *009 

30  . 74 

40 . 0 

- . 600-006 

. 2 36-006 

2.606*009 

31  .27 

10.0 

-. 100-006 

. 600-006 

3. 166*009 

39 . 40 

40 . 0 

- . 600-006 

.600-006 

2. 770*009 

34.69 

10.0 

-. 100-006 

. 100-006 

3 . 169*009 

39.67 

40.0 

-.600-006 

. 100-006 

2.849*089 

36.67 

2 0 . 0 

-.  100-006 

. 100-006 

2 . 376  *009 

29.66 

40.0 

-.  100-006 

. 1 00-086 

2 . 104*089 

26.27 

20.0 

-. 1 00-006 

.236-006 

2.668*009 

32 . 06 

40 . 0 

-.100-006 

. 2 36-006 

2 . 668*809 

3 2.07 

20 . 0 

100-006 

.600-006 

2.666*009 

33.27 

40.0 

100-006 

. 600-806 

2 . 849*089 

36.67 

20 . 0 

-.100-006 

. 100-006 

2 .691  *009 

33.69 

40.0 

-.100-006 

. 100-006 

2.932*009 

36 . 60 

20.0 

-.  2 36-006 

. 100-006 

2 . 660*009 

32 . 06 

20 . 0 

-.236-006 

.236-006 

2. 796*009 

34 . 09 

20.0 

-.236-006 

. 600-006 

2.910*009 

36.  33 

20 . 0 

-.236-006 

. 100-006 

2.941*009 

36. 72 

20.0 

-.600-006 

. 100-006 

2.666*009 

33.27 

20.0 

-.600-006 

.236-006 

2.910*009 

36.  33 

20.0 

-.600-006 

. 600-006 

3 . 0 14*009 

37. 0H 

20 . 0 

-.600-006 

. 100-006 

3 . 068*009 

38.  30 

I \BI  I Kit 


1 "’(to 

“1 

X III  " m-1 

K.  nt  i 

1 5 ■ 

ID'1  K.hI 

Imi"  1 

r-  \ 

x 

yr 

/ 1 moo  m 

* - 

**. 

At* 

.. 

Id.  0 

- . 100-00*. 

.100-40*. 

7.66(3*008 

31 .91 

20.0 

-. 100-046 

. 100-406 

1 . 44  3 *409 

60 . 40 

\ 4.0 

1 00-00*. 

.2)6-006 

1 . 066*009 

4 3.91 

2 4.0 

- . 1 00-006 

. 2 36-046 

2 .429*009 

101  .08 

id . 0 

~ . 1 00  00* 

. 604-006 

1 . 268*00-* 

62.76 

29 . 4 

-.100  046 

. 604-086 

2.699*049 

108. IS 

10.«» 

100-006 

. 100-00*- 

1 . 137*008 

66.66 

29. 4 

-.144-006 

. 100-00S 

2 .64  7*009 

110.14 

10.0 

- . 2 36-046 

. 140-006 

1 . 066*009 

43.91 

44.0 

-. 1 00-406 

. 1 04-046 

9.06  3*008 

37.72 

10.0 

- .236-006 

. 2 36-006 

1 . 691 *00  * 

19  . 36 

40 . 0 

- . 1 00-006 

. 2 36-046 

1 . 242*009 

SI  . 70 

10.0 

- . 2 36-006 

. 600-006 

2.  312*009 

96 .21 

40 . 0 

-.100  -006 

. 600-006 

1 . 489*009 

6 1 . 4S 

10.0 

-.236-006 

. 1 00-00*. 

2.  SS S*004 

106. 34 

40 . 0 

-.  1 00-40*. 

. 1 04-  80S 

1 . 669*009 

6S.  31 

10.0 

- .600-006 

. 1 00-046 

1 . 268*00  9 

62.  76 

40.4 

-. 216-006 

. 1 00-006 

1 . 242*009 

61  . 19 

10.0 

-.600-006 

. 2 36-006 

2 . 112*009 

96.2! 

40 . 0 

-.  2 36-046 

. 2 36-046 

1 .866*004 

77.21 

10.0 

-.600-006 

. 640-006 

2.662*00* 

110.  39 

40.0 

- . 2 36-046 

. 640-086 

2.067*009 

96.03 

10.0 

— .600-006 

. 100-006 

2.681*009 

111  . H 

40 . 0 

-. 2 16-446 

. 104-046 

2 . 1 30*00  * 

88. 6S 

10.0 

-.  1 00-00*. 

. 1 00-006 

1 . 317*009 

66.66 

40 . 0 

- .604-406 

. 1 04-006 

1 .489*009 

61 .96 

10.0 

-.  100-00*. 

. 2 36-006 

2.SSS*00* 

1 06. 34 

40 . 0 

- . 600-406 

. 2 36-006 

2.067*009 

86.03 

10.0 

1 00-00*. 

.600-006 

2.681*009 

111.68 

40 . 0 

- .604-006 

.600-006 

2 . 3 34  *009 

47.12 

10-0 

100-40*) 

. 100-04*. 

2. 710*009 

112.79 

40 . 0 

- . 600-406 

. 1 04-406 

2.41 4*409 

100.47 

29  . 0 

- . 1 04-006 

. 1 00-006 

8 . 462*008 

3 3.66 

40.0 

-. 104-006 

. 1 04-046 

1.669*009 

66.  31 

.00 

-. 1 00-006 

. 2 J6-006 

1 . ] 08*009 

46.12 

40 . 0 

-. 100-006 

. 2 36  -046 

2 . 1 30*004 

88 . 66 

.0  0 

100-006 

.600-006 

1 . 3 30*009 

66.  3 7 

40 . 0 

-. 100-006 

. 600-806 

2.414*009 

100.47 

29 . 0 

-. 100-006 

. 100-00*. 

1 .403*009 

68 . 40 

40 . 0 

-. 140-406 

. 100-046 

2.601 *009 

104.07 

29 . 0 

- . 2 36-006 

. 1 00-006 

1 . 108*009 

46.12 

00  0 

-.  2 36-006 

. 2 36-006 

1 . 772*009 

7 3.73 

.0.0 

-. 2 16-006 

.*>00-006 

2 . 388*009 

99 . 40 

29  0 

- . l 36-006 

. 160-40*. 

2 .429*009 

101.08 

29.9 

-.600-006 

. 1 00-006 

1 . 3 30*009 

66.37 

29 . 9 

- . 644-006 

. 2 36-006 

2.  388*009 

99 .40 

.0.0 

-.600-006 

. 600-006 

2.SS2*009 

106.22 

<4.0 

-.604-006 

. 100-00*. 

2.699*009 

108.16 

4 


I A HI  I Km  i 


i.  = 7ih)  a|  ^ :k»  x nr'1  m_l 


Hi' 

it))  = 15  X 

lO"1  Kid 

1m  m i 

= 15  X 

I0~3  R.id 

/ 

, = 1500 

111 

♦»/ 

*. 

01  1 * 

lH/.t  M/ 

ATNU 

If  0 

-.  100-006 

. 100-006 

1 .660*009 

2 . 00 

20 . 0 

- 

100-006 

. 100-006 

1 .614*009 

2.02 

11.0 

100-006 

. 2 16-006 

1 . 69  7 *00* 

2.12 

20 . 0 

- 

1 00  -00  S 

.2)6-006 

1 .694*009 

2.11 

10.0 

-.  1 00-006 

. 600-006 

1 . 710*009 

2.14 

20.0 

- 

100-006 

. 600-006 

1 .7)2*009 

2.16 

10  0 

-. 1 00-00  6 

. 100  -00' 

1 . 714*004 

2.14 

20 . 0 

- 

100-006 

. 100-006 

1 . 742*009 

2.  8 

10.0 

-.2  )6-#06 

. 100-006 

1 .647*009 

2.12 

40.0 

- 

100-006 

. 1 00-006 

1 . 1 77*004 

1 .47 

10.0 

-.2  36-006 

. 2 36-006 

1 . 727*009 

2.16 

40 . 0 

- 

1 00-006 

. 2 36-006 

1 . 303*009 

1.63 

10  0 

- . 236-006 

.600-006 

1.741 *009 

2.17 

40 . 0 

- 

110-086 

.600-806 

1 . 369*009 

1 .71 

10.0 

236-006 

. 100-006 

1 . 744*009 

2.18 

40.0 

- 

100-006 

.100-006 

1 . 387*009 

1.73 

10.0 

- . 60? -006 

. 1 00-006 

1 . 710*009 

2.14 

40 . 0 

- 

2)6-006 

. 100-806 

1 . 30  3*089 

1.6  3 

10.0 

-. b00-00b 

. 2 36-006 

1.741 *009 

2.17 

40 . 0 

- 

2)6-006 

.2)6-886 

1 .460*009 

l .02 

10.0 

- . b00-00b 

. 600  -006 

1 . 766*009 

2.19 

40 . 0 

- 

236-806 

.600-006 

1 .542*009 

1.93 

10.0 

-. 600-006 

. 100-006 

1 . 760*004 

2.20 

40.0 

- 

236-006 

. 100-006 

1 . 666*009 

1 . 96 

10.0 

100-006 

. 1 00-006 

1 .714*009 

2.14 

40.0 

- 

600-006 

. 100-006 

1 . 369*009 

1 . 71 

10.0 

100-006 

. 2 36-006 

1 . 744*009 

2.18 

40 . 0 

- 

600-00 b 

.2)6-006 

1 .642  *089 

1.93 

10.0 

-.  100-006 

. 600-006 

1 . 768*009 

2.20 

40 . 0 

- 

600-006 

.680-006 

1 .6  36*009 

2.04 

10.0 

1 00-00*) 

. 100-006 

1 . 762*009 

2.20 

40 . 0 

- 

600-006 

. 100-006 

1 .661  *009 

2.07 

2 0 . 0 

-. 1 00-006 

. 100-006 

1 . 604  *009 

1 . 88 

40 . 0 

- 

100-006 

. 100-006 

1 . 387*009 

1.73 

20  . 0 

-. 100-00b 

. 2 36-006 

1 . 67  3*009 

1.96 

40.0 

- 

100-006 

. 2 36-086 

l . 666*009 

1 .96 

20.0 

1 00 -00  b 

. 600-006 

1 .606*009 

2.00 

40.0 

- 

100-006 

. 600-006 

1.661 *009 

2.0  7 

20 . 0 

-.100-0*6 

. 100-006 

1 .614*009 

2.02 

40.0 

- 

100-006 

. 100-006 

1 .688*009 

2.11 

20.0 

-.236-006 

. 100-006 

1 .673*009 

1 .96 

20.0 

-.236-006 

. 2 36-006 

1 . 648*009 

2.06 

20.0 

-.2)6-006 

.600-006 

1 .604*009 

2.10 

20.0 

- . 2 36-006 

. 100-006 

1 .694*009 

2.11 

20.0 

— .600-006 

. 100-006 

1 . 606*009 

2.00 

20.0 

• . 600-006 

.236-006 

1 .604*009 

2.10 

20.0 

-.600-006 

.600-006 

1 . 722*009 

2.16 

2 0.0 

-.600-006 

. 100-006 

1 . 732*009 

2.16 

43 


I \H I I '(ni 


1 

700 

(‘l 

2} f>  X 1 0 in 

1 

Ko 

(Hi  1 5 > 

nr1  k.hI 

Imdl  1 

50 

x II)-'  K.ul 

'\ 

1 500 

m 

10.0 

- . 100-006 

. 1 00-00* 

1 .6*4*004 

6.41 

20.0 

-.  100-00* 

. 1 00-006 

.611 *009 

6.  71 

10.0 

. 103-00* 

. 2 16-00* 

1 .644*004 

7 .0* 

20 . 0 

-.  1 00-00* 

. 2 16-006 

.642*004 

7.04 

10.0 

- . 1 00-00* 

. *00 -00* 

1 . 70  7 *004 

7.11 

20.0 

100-00* 

. *00-00* 

1.710*009 

7.20 

10.0 

-. 1 00-00* 

. 100-00* 

1.71 1*004 

7.12 

20.0 

-.  100-00* 

. 100-00* 

i.740*009 

7.24 

10.0 

- . 2 lb-00* 

. 1 0 0 - 0 0 h 

1 .644*004 

7.0* 

40.0 

-.  100-00* 

. 100-006 

.174*004 

4.89 

10.0 

- . 2 16-00* 

. 2 16-00* 

1 . 72**004 

7.  1H 

40 . 0 

100-006 

. 2 16-00* 

. 1 0 1 ♦ 0 0 9 

.41 

10.0 

-.2  16-00* 

. * 0 0 -0  0 * 

1 . 7 14*004 

7.24 

40 . 0 

-.  100-006 

.600-006 

1 . 166*009 

S .64 

10.0 

- . 2 16-0  0b 

. 100-00* 

1 . 742*004 

7.2* 

40.0 

100-006 

. 1 00-00* 

i . 184*009 

*.76 

10.0 

- ' 00-00* 

. 1 00-00* 

1 . 707*004 

7.11 

40.0 

-.2  16-006 

. 100-00* 

1 . 101*009 

*.41 

10.0 

-.600-00* 

. 2 16-00* 

1 . 7 14*004 

7.24 

40 . 0 

-.216-006 

. 2 16-00b 

i . 4*7*009 

6.0  7 

10.0 

-.  *00-006 

.600-00* 

1 . 7*1*004 

7.24 

40.0 

- . 2 16-006 

.600-00* 

1 . *40*009 

6.41 

10.0 

- . *00-00* 

. 100-00* 

1 . 7*7*004 

7.11 

40 . 0 

- . 2 16-006 

. 100-00* 

1 .*61*009 

6 . * 1 

10.0 

100-00* 

. 100-00* 

1 .711*004 

7.12 

40 . 0 

- .600-006 

. 100-00* 

l . 166*00  4 

*.69 

10.0 

- . 1 00-00* 

. 2 16-00* 

1 . 742*004 

7.2* 

40.0 

-.600-006 

. 2 16-006 

1 . *40*004 

6.41 

10.0 

-.100-00* 

. *00-00* 

1 . 7*7*004 

7.  11 

40.0 

- .600-006 

. 600-00* 

1.6)  1*009 

6 . 8? 

10.0 

100-00* 

. 100-00* 

1 . 7*0*004 

7.11 

40 . 0 

- . 600-006 

. 100-00* 

. 6*9*009 

6 . 90 

00.0 

1 00-00* 

. 1 0 0 - 0 0 b 

1 . *01 *004 

6.24 

40.0 

-.  100-00* 

. 100-006 

i . 184*004 

*.  76 

20.0 

100-00* 

. 2 16-00* 

1 . *70*004 

6.  *1 

40.0 

100-00* 

. 2 16-006 

. *6  1*009 

6 . * 1 

20.  0 

1 00-006 

. 600-00* 

1 . 60  1*009 

6.67 

40 . 0 

100-00* 

. 600—006 

.**4*00  * 

* . 40 

00.0 

-. 100-00* 

. 1 00-00* 

1.611 *004 

6.  71 

4 0.0 

-.100-00* 

. 100-00* 

.686*004 

7 .02 

. 0 . 0 

- . 2 16-00* 

. 100-006 

1 . *70*004 

6.  *3 

<’0.0 

- . 2 16-00* 

. 2 16-0  06 

1 .646*004 

6 . 8* 

20 . 0 

-.2  16-00* 

.600-00* 

1 . 602*004 

7.00 

.0.0 

-.2  16-00* 

. 100-00* 

1 . 642*004 

7.04 

2 0 . 0 

-.*00-00* 

. 100-00* 

1 .60  1*004 

6.67 

20 . 0 

-.*00-00* 

. 2 lb-006 

1 .682*004 

7.00 

20.0 

-.*00-00* 

.600-00* 

1 .720*004 

7.16 

20 . 0 

-.*00-00* 

100-00* 

1 . 7 10*004 

7.20 

44 


I AH  I I ,*(o) 


I = 700.  a | = ,23ft  X 10-<1  m_1 


Rt 

< 0 ) = 1 5 X 

10'1  R.ul 

lm<«  l 

II 

X 

I0~:  Rail 

/ 

',|  = 1500. 

111 

„ 

f MA  * Mi 

AT 6 . 

IB  . 

Of  l M 

- 1 

*3 

f Mf<(  M / i 

i'iB/v  - 

19.9 

- . 1 00-906 

. 100-006 

1 .6  16*009 

20.43 

20.0 

-.  100-006 

. 100-006 

1 . 688*004 

19.8) 

ii.i 

-.  100-006 

.236-00' 

1 .669*004 

2 0.84 

20.0 

-.  1 00-006 

. 2 36-006 

1 .672*009 

20.88 

.9  i 

-.  100-006 

. 600-006 

1 .664*009 

21.0  3 

20.0 

-.  1 00-006 

. 600-006 

1 .713*009 

21  . 38 

»#.# 

-.  199-996 

. 100-006 

1 .660*009 

21.08 

20 . 0 

-.  100-006 

. 1 00-006 

1 . 724*009 

21.62 

19.0 

-.236-996 

. 100-006 

1 .669*009 

20.64 

40.0 

-.  1 00-006 

. i 00-006 

1.163*004 

14.40 

10.0 

-.236-996 

.236-006 

1 . 704  *004 

21  .27 

40 . 0 

-.  100-006 

. 2 36-006 

1 . 280*009 

16.98 

10  .i 

-.2 3b-90b 

.600-006 

1 .719*004 

21.47 

40.0 

-.  100-006 

.600-006 

1 . 346*004 

16.80 

10.0 

-.2 36-906 

. 100-006 

1 . 724*009 

21.62 

40 . 0 

-.  100-006 

. 100-006 

1 . 364*004 

17.03 

10.0 

-.600-996 

. 1 00-006 

1 .664*009 

21.0) 

40 . 0 

- . 2 36-008 

. 100-006 

1 . 280*004 

16.98 

10.0 

-.600-906 

. 2 36-006 

1 . 719*009 

21.47 

40 . 0 

-.236-006 

. 2 36-006 

1 . 438*009 

17.96 

10  0 

-.600-996 

. 600-006 

1 . 7 36*004 

21.67 

40 . 0 

-.  2 36-006 

.600-006 

1 .622*004 

19.00 

10.0 

- .600-906 

. 1 00-006 

1 . 740*004 

21  . 72 

40 . 0 

-.  236-006 

. 1 00-006 

1 . 646*009 

19.29 

10.0 

-.  199-006 

. 1 00-006 

1 .688*004 

21  .06 

48 . 0 

- .600-006 

. 1 00-006 

. 346  *00  4 

16.80 

<9.0 

-.  199-906 

.236-006 

1 . 724*009 

21  . 62 

40 . 0 

-.600-006 

. 2 36-006 

1 . 622*004 

19.00 

19.9 

-.  190-906 

.600-006 

1.740*009 

21.  72 

40 . 0 

- . 60  0-006 

.600-006 

1 .616*009 

24.17 

19.9 

-. 1 99-906 

. 100-006 

1 . 744*004 

21.77 

40 . 0 

- .600-006 

. 1 00-006 

1 .642*004 

20.60 

20.  i 

-.  100-906 

. 100-006 

1 . 4 7 3*009 

18.39 

40.0 

-.  1 00-006 

. 100-006 

1 . 364  *009 

1 7.0  3 

20 . 9 

-.  1)30-006 

. 236-006 

1 . 646*004 

19.29 

40.0 

- . 100-006 

.236-906 

1 . 646*009 

19. 29 

29 . 9 

100-006 

. 600-006 

1 . 679*004 

19.72 

40 . 0 

-.  1 00-006 

. 600-006 

1 .642*009 

20 . 60 

29. 9 

100-006 

.100-006 

1 .689*004 

19.83 

40.0 

100-006 

. 100-006 

1 .664*009 

20  . 84 

29.9 

-.236-006 

. 100-006 

1 . 646*004 

19.29 

29.9 

-.236-006 

.236-006 

1 . 624*004 

20.28 

29.9 

-.236-006 

.600-006 

1 .662*009 

20.76 

29.0 

-.236-006 

. 100-006 

1 .672*004 

20. 68 

20.9 

-.600-006 

. 100-006 

1 . 679*004 

19.72 

29.9 

-.600-006 

.236-006 

1 .662*004 

20.76 

20.9 

-.600-006 

. 600-006 

1 . 702*004 

21.26 

20 . 9 

-.600-006 

. 100-006 

1.71  3*009 

21.38 

45 


I AH  I I ( p i 


i = 7ou  o , :o.  x io  m"' 


Kc 

uM  ; 1 5 x 

nr1  k id 

ImifJ ) 

.50  / 

10"  - Rad 

/,  1500. 

111 

. .*  «* 

- . i 83-00* 

. 100-006 

7 . 440*00-3 

10 . 96 

20.0 

-.100-306 

. 103-006 

1 . 329*009 

16/.  ■ 1 

66.  11 

10.0 

- . 1 00-03b 

. 2 16-006 

1 .312*009 

42.11 

20 . 0 

-. 100-006 

. 2 36-006 

1 . 4 72*009 

61.26 

10.  * 

- . t 33-30*. 

. 600-006 

1 .207*009 

60.21 

20.0 

-.100-006 

. 600-006 

1 .633*009 

6 3.78 

1#.** 

100-006 

. 1 00-00  S 

1.270*009 

62.84 

20. 0 

-.100-006 

. 1 00-006 

1 .649*009 

64.47 

10.3 

-.2  J6-036 

. 100-006 

1 .012*009 

42.11 

40.0 

-. 103-006 

. 100-006 

8. 746*038 

36.40 

10.3 

- . 2 36-386 

.236-006 

1 .460*009 

61.60 

40 . 0 

-. 100-006 

. 2 16-006 

1 .076*009 

44.73 

: *.  3 

- . 2 16-006 

.600-006 

1 »SlS+tV9 

63.04 

40 . 0 

-. 100-006 

.600-006 

1 . 143*009 

47  ,6 

l «.** 

-.2  16-006 

. 100-00S 

1 .624*009 

6 1.41 

40.0 

-. 100-006 

. 1 00-006 

1 .162*009 

4 H . 34 

10.3 

-.600-006 

. 100-006 

l . 20  7*009 

60.21 

40. 0 

- . 2 36-00b 

. 100-006 

1 .076*009 

44.73 

it.* 

- . 600-006 

. 2 16-006 

1 . 616*009 

61.04 

40 . 0 

-.  2 36-006 

. 2 36-0 0b 

1 .240*009 

61.69 

10.0 

- . 600-006 

. 600-006 

1 .661*009 

h4 . 66 

40 . 0 

-.  2 36-006 

.600-006 

I . 331*009 

66.  38 

10.0 

-.600-006 

. 100-00S 

1.661 *009 

64 . 9b 

40.0 

-.236-006 

. 1 00-006 

1 . 367*009 

66.46 

i 0 . 

- . 100-006 

. 100-006 

l . 270*009 

6 2.84 

40.0 

-.600-006 

. 100-006 

1 . 14  1*009 

47.66 

19.* 

-. 108-336 

. 2 36 -006 

1 . 624*009 

63. 4 J 

40.0 

-.600-006 

.236-006 

1 . Ill *009 

66.  18 

ie.  * 

103-036 

. 600-006 

1 .661*009 

b 4 . 96 

40.0 

-.600-006 

. 600-006 

1 .436*009 

69. 78 

10.0 

100-006 

. 100-3*6 

1.671*009 

66.  37 

40 . 0 

-.600-006 

. 100-006 

1 . 467*009 

61.01 

20. 0 

- . 1 00-006 

. 100-006 

7.011 *00  8 

32.61 

40 . 0 

-.100-006 

. 100-006 

1 . 162*009 

49  . 14 

19.9 

-.  100-006 

. 2 16-0  06 

1 .061*009 

4 4.  1 6 

40.0 

-.100-006 

.236-006 

1 . 367*009 

66.46 

2 0.3 

-. 100-006 

.600-006 

1 . 263*009 

62.67 

40 . 0 

-.100-006 

. 600-006 

1 . 467*009 

61.01 

.0  0 

-.  1 00-006 

. 100*006 

1 . 329*009 

66. 31 

40 . 0 

-.100-006 

. 100-006 

1 . 498*009 

62.  14 

20.0 

-.2  16-006 

. 100-006 

1.061 *009 

44.16 

.0.0 

-.2  16-006 

. 2 16-006 

1 .402*009 

68.36 

00 . 0 

-.2  16-006 

. 600-006 

1 .467*009 

60.64 

20.0 

-.2  16-006 

. 1 00-336 

J . 4 72  *009 

6 1.26 

Si. 9 

600-006 

. 1 00-006 

1 .26  1*009 

62.67 

20.1 

- .600-006 

. 2 16-006 

1 . 46  7*009 

60.64 

00 . 0 

- . 600-006 

. 600-006 

1 .616*009 

63.11 

00 . 0 

- . 600-006 

. 1 00-006 

1.631  *009 

63.78 

4<> 


7 CONCLUSIONS 


lo  use  the  RIM  as  a search  scheme  lor  modal  solutions,  the  original  modal  equation 
(equation  3)  expressed  in  terms  ol  plane  wave  reflection  coefficients  R[)<0)  and  R[)(tf)  was 
modified  m such  a was  as  to  ef  fectively  move  its  poles  to  the  point  infinity  in  the  complex 
plane  One  reason  lor  this  modification  is  clearly  depicted  in  figures  I (land  II.  figure  10 
exhibits  "phase  line  extinction"  associated  with  poles  which  occur  within  the  search 
rectangle  "Phase  line  extinction”  makes  it  possible  lor  the  RIM  to  miss  solutions  within 
the  search  rectangle  Figure  I I shows  the  much  relaxed  phase  contour  structure  which 
results  with  all  poles  moved  to  the  point  infinity.  For  such  phase  contour  behavior,  the 
R!  \!  will  pick  up  all  modes  within  the  search  rectangle.  I he  modification  ot  the  original 
modal  equation  has  one  serious  disadvantage  in  the  sense  that  numerically  large  factors 
occurring  in  the  numerators  and  demoninators  of  quantities  which  occur  in  the  original 
mode  equation  and  which  can  be  simply  canceled  out  in  the  ratio,  occur  as  products, 
rather  than  ratios,  in  the  modified  modal  equation.  I his  disadvantage  can  result  in 
overflow  problems  in  the  microwave  range. 

For  the  case  study  ot  the  trilinear  layer  ol  figure  R.  implementation  of  a single 
precision  version  ol  the  RI  M on  a Univac  1110  was  straightforward  and  successful  for  the 
to-.  I TF.  and  520-MHz  cases  examined.  Figure  (>  is  a representation  ol  the  RIM  applied  to 
the  52U-MII/  case,  this  being  an  interesting  example  ol  this  lower  frequency  domain.  Figure 
')  illustrates  the  interesting  feature  that  the  RIM  does  not  necessarily  locate  roots  sequen- 
tially according  to  the  magnitude  ot  Re(0  ). 

Overflow  was  a problem  for  the  3.3-011/  case,  although  it  would  not  have  been  had  a 
double  precision  version  ol  the  RI  M been  used.  Function  switching  (see  section  5 I and 
search  rectangle  subdivision  (see  section  6)  have  been  utilized  to  overcome  the  overflow 
problems  in  this  case  m such  a way  that  computer  run  times  were  comparable  with  those  of 
the  more  conventional  Newton-Raphson  scheme.  Unlike  the  latter,  of  course,  the  RIM 
assures  the  location  ol  all  zeros  (or  modes)  within  the  search  rectangle 

Unfortunately . f unction  switching  along  with  search  rectangle  subdivision  is  very 
difficult  to  automate  (this  would  be  especially  true  for  nuiltisegmented  profiles).  I hus.  the 
RIM  appears  to  have  its  greatest  utility  tor  tropospheric  propagation  when  function  switch- 
ing and  search  rectangle  subdivision  can  be  altogether  avoided.  For  this  reason  a set  of  tables 
(tables  3(a)  through  3(p)>  has  been  prepared  as  a guide  to  the  waveguide  program  user.  The 
tables  show  an  estimated  maximum  frequency.  1MAX.  below  which  function  switching 
should  not  occur  lor  several  trilinear  layer  models  and  various  search  rectangles.  I he  method 
ol  calculation  should  assure  a rather  conservative  estimate  of  1 MAX  lie.  switching  should  not 
occur  until  frequencies  somewhat  higher  than  (MAX  are  encountered).  Since  there  are 
many  practical  instances  when  FMAX  is  on  the  order  of  5 (ill/  or  more,  we  recommend 
implementing  the  RI  M with  a nuiltisegmented  waveguide  model 
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APPENDIX  A DERIVATION  Ol  THE  POLE-1  REE  C ONTINUOUS 
VERSION  Ol  THE  MODAL  EQUATION 

A pole-free  discontinuous  version  of  the  modal  equation  is  given  by  equations  (24 1 
through  (27)  of  section  5 In  this  appendix  we  sketch  the  derivation  ol  the  pole-tree  contin- 
uous version  given  by  equations  ( 28 ) through (41  ) of  section  V The  essence  of  the  deriva- 
tion depends  upon  the  asy  mptolic  behavior  of  the  denominators  D/2.  I ) A 2 . and  DBI  given 
In  equations  (4 ).  (5 ).  and  (17)  of  section  5 since  the  switching  Irom  one  function  to  anoth- 
er occurs  only  it  the  unaginai y part  ot  one  or  all  ol  the  quantities  p p I ).  pil  2 ).  and  p | ( I i is 
sufficiently  large  (as  determined  by  L in  the  inequality  ( 2 3 1 ( . E'urthermore.  we  also  know 
from  the  physics  of  the  problem  that  I in  lpp  I ).  pa(  2).  p j < I ))  > 0 It  is  assumed  that  the 
switch  front  one  function  to  the  other  occurs  only  when  the  modified  refractive  index  is 
sufficiently  negative.  Mathematically  this  condition  is  formulated  by  requiring  that  the 
relevant  quantity  or  quantities,  Arg  ( P3<  I )).  Arg  ( ps(  2 )),  Arg  ( p | ( I ))  be  greater  than  2 3 rr . 

I hus.  summarizing  to  this  point,  we  are  concerned  with  the  asymptotic  development  of  the 
denominators  I)/2.  I)A2.  and  1)B1  in  the  region  of  the  complex  plane  bounded  by  the  radi- 
al 2 4 rr  and  the  negative  real  axis.  The  appropriate  asymptotic  expansions  ol  the  modified 
Hankel  function  of  order  one-third  are  therefore: 


-14  t S \ 

hj  tql  ~ 0q  exp  j-^iq-  “ - p jti  J 

• 14  / ^ t ■>  S \ 

h,  1 1 ) ~ i 0q  exp  i <p  --pm) 

- 1 4 / 1 } s 11  \ 

hs  <qt  " 0 4 exp  I ~ i q-  " + p m j 

1 14  / ^ ^ *>  1 1 \ 

lis  Iql  " I0‘1  exp  { ~ i q-  -+  — m ) 


(A  I) 

( A 2 ) 

( A 4 I 


(A  4) 


I rom  these  formulas,  the  asymptotic  developments  of  the  denominators  D/2,  DA2. 
and  DBI  given  by  equations  (4),  (5).  and  ( 17)  of  section  5 are: 

D/2  - id<P|  (4))1  4 exp  i (p|(3)^  2 * + p m ) (A. 5) 

DA 2 - 20(P2<2))-*  4 exp  i ( pM  2 ))■'  - - p m ) ( A (>) 

DBI  ' P ( p j ( I If1  4 ( I - Vq  ( s-  - «|  II ) 1 ~j  exp  (-=-i<P|<  I ))4,  “ - p rri)  ( A.  7) 

Now  observe  that  I | given  by  equation  (24)  of  section  5,  differs  from  I 4 given  by 
equation  1 27 ) by  the  laetors  DZ2  X DA2  X 1)7.0  X DB  I I hus.  one  way  they  can  be  made 
continuous  at  the  switching  boundary  is  to  multiply  I | by 

( DZ2 )~ 1 X ll)A2f 1 X (DBI  )_1  (A.8) 

where  the  denominators  are  the  asymptotic  expressions  given  m equations  ( A. 5 ) through 
( A 7)  and  to  multiply  I 4 by  1)7,0.  Next  E s as  given  by  Equation  ( 25 ) ol  section  5 can  be 
made  continuous  at  the  switching  boundary  with  this  newly  formed  I 4 by  multiplying  the 
I s by  ( I)B  I r * Similarly.  E } given  by  equation  ( 2(>  I of  section  5 can  be  made  continuous 


4‘) 


with  the  newly  formed  I 4 .it  the  switching  boundary  hy  multiplymj!  I 3 hy  <l)/2l  1 X 
1 1)  \ 1 r'  1 inalK  . the  given  in  eci  nations  ( 2X ) through  (3  I I ol  sect  ion  5 are  obtained 
In  multiply  me  the  new  set  ol  1,'s  through  by  the  laelor.  (p|  ( I tr1  '4ll-.()lS- 
- a | 1 1 1 1 - 1 
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